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Abstract 

This  final  report  and  draft  guidance  document  focuses  on  assessing  the  level  of  site 
characterization  needed  to  support  quality  decisions  regarding  remedial  strategies  and  long  term 
stewardship  of  contaminated  sites  using  a  flux  and  mass  balance  based  approach.  Site 
characterization  efforts  are  aimed  at  understanding  the  link  between  dense  nonaqueous  phase 
liquid  (DNAPL)  source  zones  mass  discharge  under  natural  hydrological  conditions  and 
modified  as  a  result  of  source  zone  treatment.  The  transition  between  forward  and  back- 
diffusion  was  also  explored.  This  document  provides  an  overview  of  methodologies  that  can  be 
used  to  characterize  source  zone  mass  discharge  functions  and  how  this  controls  plume  response 
and  back  diffusion. 

Prediction  of  the  relationship  between  DNAPL  source  zone  changes  and  down-gradient  plume 
response  is  critical  to  making  informed  site  management  decisions,  especially  those  related  to 
remedial  actions.  Thus,  it  is  vital  that  source  zone  and  plume  characterization  be  conducted 
within  a  framework  that  is  consistent  with  appropriate  predictive  models.  The  project  focused  on 
demonstration  of  effective  field-scale  approaches  that  provide  vital  links  between 
characterization,  prediction,  and  decision  making. 

Specific  objectives  of  this  project  were  to:  1)  Develop  source-strength  functions  for  site 
management  purposes  using  existing  historical  site  data  supplemented  with  limited  flux-  and 
core-based  sampling;  2)  Extend  to  the  field-scale  our  ability  to  predict  DNAPL  source  depletion 
through  dissolution,  based  on  a  priori  characterization  of  the  source  zone  architecture;  3) 
Characterize  the  near-source  plume  response  to  source-mass  depletion  at  selected  field  sites  in 
order  to  provide  the  understanding  needed  to  predict  long-term  plume  responses;  4)  Link  the 
characterization  of  the  near-source,  short-term  responses  to  likely  long-term  behavior  of  the 
dissolved  plume  including  back  diffusion;  and  5)  Provide  recommended  guidance  on  the  level  of 
source  zone  characterization  needed  to  adequately  predict  source- strength  functions  and  plume 
response. 

The  research  effort  developed  source- strength  functions,  Mo(t),  using  existing  historical  data  at 
field  sites.  We  selected  five  primary  field  sites  and  up  to  ten  secondary  sites  with  available  sets 
of  high  quality  historical  data.  At  the  primary  sites,  additional  data  was  collected. 
Characterization  of  the  source-strength  function  were  based  on  a  quality  historical  record  of 
concentration  and  head  measurements,  particularly  in  the  near-source  region.  Lor  each  site 
considered,  currently  available  data  alone  was  used  to  provide  the  best  possible  description  of  the 
source  function.  Lollowing  completion  of  this  exercise  on  secondary  sites,  primary  sites  were 
selected.  Sites  were  evaluated  to  identify  what  type  of  additional  data  will  be  obtained  including 
flux-based  measurements  such  as  integral  pump  tests  and  passive  flux  meter  (PPM)  well  screen 
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profiles.  The  data  collected  from  the  primary  sites  was  used  to  improve  the  estimated  source- 
strength  function  and  associated  error  estimates. 

Advanced  characterization  of  the  DNAPL  architecture  included  current  source  zone  mass 
(Mo),  contaminant  fluxes,  J(x,y,z),  and  mass  discharge  rate,  Md;  high-resolution  DNAPL 
saturation  characterization,  SN(x,y,z);  and  trajectory-integrated  DNAPL  saturation.  Mass  flux 
measurements  made  at  down-gradient  control  planes,  and  tracer-based  trajectory-integrated  data 
are  Lagrangian  characteristics  that  contain  information  representing  the  actual  flow  paths  taken 
by  water  traversing  the  source  zone.  Point-based  spatial  data  such  as  soil  cores  contain  different 
information,  and  both  data  types  will  be  used  here  to  predict  dissolution. 

At  the  Alameda  California  primary  site  selected  for  further  evaluation,  the  down-gradient 
response  to  changing  source  zone  strength  was  quantified  through  a  controlled  field  experiment. 
The  experiment  was  conducted  by  containing  a  portion  of  an  existing  source  zone  using  three 
injection  wells  and  four  extraction  wells  forming  a  hydraulic  control  cell.  Within  this  flow  cell, 
the  source- strength  function  behavior  under  enhanced  groundwater  flow  conditions  was 
observed.  Reactive  and  nonreactive  tracers  were  used  to  characterize  transport  conditions  within 
the  flow  cell  including  diffusion  into  lower  conductivity  zones  with  the  aquifer.  The  subsequent 
response  of  the  contaminant  plume  to  a  source- strength  function  perturbation  was  monitored  for 
“back  diffusion”  from  lower  conductivity  zones  or  DNAPL  contaminated  zones. 

The  field  data  collected  in  this  project  was  used  to  evaluate  several  approaches  for  estimating  the 
source  strength  function,  and  the  uncertainties  associated  with  the  estimated  model  parameters. 
Cost-effective  approaches  to  DNAPL  site  characterization  allows  Department  of  Defense  (DoD) 
users  and  site  managers  to  more  accurately  assess  the  benefits  of  costly  aggressive  source  zone 
treatment  technologies.  Accurate  assessment  of  near-source  flux  changes  and  the  magnitude  of 
"back  diffusion"  expected  will  help  DoD  by  avoiding  costly  remedial  efforts  with  inadequate 
benefits. 

The  results  of  the  project  can  be  used  by  DoD  site  managers  and  consultants  to  better 
characterize  sites  in  order  to  make  informed  decisions  regarding  implementing  aggressive  source 
zone  remedial  efforts.  Guidelines  on  site  characterization  efforts  aimed  at  understanding  the  link 
between  the  site  DNAPL  source  strength  function  under  natural  hydrological  conditions  and  as 
modified  by  source  treatment.  These  guidelines  provide  an  established  methodology  for  defining 
the  source  strength  function  and  how  this  function  defines  predicted  plume  responses.  The  ability 
to  better  characterize  these  relationships  can  reduce  site  management  and  clean-up  costs  by 
enabling  quantitative  prediction  of  the  benefits  of  costly  source  zone  treatment.  In  addition,  this 
project  developed  and  disseminated,  through  journal  articles,  recommended  guidelines  and 
workshops,  an  approach  for  quantifying  "back  diffusion"  processes  at  contaminated  sites. 
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1.  Objectives 

Specific  objectives  of  this  project  were  to:  1)  Develop  source-strength  functions  for  site 
management  purposes  using  existing  historical  site  data  supplemented  with  limited  flux-  and 
core-based  sampling;  2)  Extend  to  the  field-scale  our  ability  to  predict  DNAPL  source  depletion 
through  dissolution,  based  on  a  priori  characterization  of  the  source  zone  architecture;  3) 
Characterize  the  near-source  plume  response  to  source-mass  depletion  at  selected  field  sites  and 
within  an  isolated  section  of  the  aquifer,  to  provide  the  understanding  needed  to  predict  long¬ 
term  plume  responses  including  back  diffusion;  4)  Link  the  characterization  of  the  near-source, 
short-term  responses  to  likely  long-term  behavior  of  the  dissolved  plume;  and  5)  Provide 
recommended  guidance  on  the  level  of  source  zone  characterization  needed  to  adequately  predict 
source-strength  functions  and  plume  response. 

2.  Background 

a.  Development  of  Mass  Flux  Based  Evaluation  of  DNAPL  Contaminated 
Sites 

The  characterization  of  DNAPL  contaminated  sites  is  shifting  the  focus  to  techniques  that  can 
lead  to  a  flux  based  assessment  of  both  contaminants  and  groundwater,  rather  than  focusing 
primarily  on  concentration  measurements.  Llux  data  provides  the  ability  to  integrate  information 
leading  to  mass  discharge,  source  mass  balances,  and  changes  in  mass  balance  over  time.  The 
primary  challenge  with  DNAPL  sites  is  generated  by  the  density  contrast  between  the  immiscible 
phase  and  water.  DNAPL  will  migrate  below  the  water  table,  and  become  trapped  on  media 
interfaces.  The  migration  of  DNAPL  results  in  varying  distributions  within  the  subsurface,  and 
challenges  our  ability  to  implement  remedial  methods. 

Site  characterization  must  begin  with  site  assessment,  addressing  the  type  of  contaminants,  the 
concentration  of  contaminants  in  water  samples,  and  the  spatial  extent  of  the  contaminated  area. 
Lor  each  contaminant  there  is  a  regulatory  maximum  contaminant  level  (MCL).  If  a  NAPL  is 
suspected,  the  concentrations  may  be  compared  to  the  solubility  limit  of  the  NAPL  mixture 
present  in  the  aquifer.  The  extent  of  the  source  zone  is  then  determined,  and  the  site  is  divided 
into  a  source  zone  and  a  plume.  Remedial  technologies  are  assessed  for  either  source  zone  or 
plume,  depending  on  the  further  characterization  of  the  zones.  Remediation  is  then  optimized  for 
cost  effectiveness  and  time  effectiveness. 

The  source  zone  and  the  plume  are  divided  distinctly  within  a  model;  however  there  may  be 
secondary  source  zones  as  well  as  back-diffusion  of  mass  that  may  extend  the  plume. 
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Figure  1.1  Conceptual  site  with  source  zone  and  plume. 


Contaminant  mass  flux  and  integrated  mass  discharge  have  been  defined  and  characterized  for 
sites.  Contaminant  mass  flux  has  units  of  mass/area/time,  and  is  the  mass  of  contaminant  moving 
across  a  unit  area  of  porous  media  perpendicular  to  groundwater  flow.  Integrated  mass  flux  is  the 
total  mass  moving  across  an  area,  over  time.  A  control  plane  may  be  defined  as  the  cross  section 
the  total  contaminant  mass  is  crossing,  and  generally  covers  the  entire  plume.  The  units  of 
measurement  for  integrated  mass  flux  are  mass/time.  Characterizing  mass  flux  and  mass 
discharge  becomes  useful  while  investigating  attenuation  processes,  whether  natural  or 
enhanced,  remedial  assessments,  measuring  loads  to  receiving  bodies,  and  risk  assessment. 

Spatial  and  temporal  methods  of  measurement  can  be  used  to  quantify  mass  flux  and  integrated 
mass  discharge.  Point-scale  techniques  use  multilevel  samplers  and  passive  flux  meters  within 
screened  wells,  and  both  methods  help  to  define  mass  flux.  The  mass  flux  data  can  then  be 
integrated  to  obtain  mass  discharge.  Integral  pump  tests  can  also  be  implemented  for  large 
volumes  of  water  being  pumped  from  the  aquifer  from  fully  screened  wells.  Integral  pump  tests 
lack  spatial  scales  of  measurement,  but  offer  a  “plume  core.”  Each  method  is  investigated  to 
determine  the  usefulness  for  site  characterization. 

Flux  measurements  have  been  developed  in  order  to  better  assess  contaminated  DNAPL  sites. 
The  primary  advancement  has  focused  on  the  relationship  between  the  contaminant  mass  present 
in  the  NAPL  source  zone,  and  the  contaminant  mass  flux  (source  strength)  which  forms  the 
plume.  Site  management  strategies  are  now  using  the  fundamental  knowledge  based  on  mass 
flux,  mass  discharge,  and  mass  balance.  Site  age,  defined  as  the  relative  mass  remaining  in  the 
source  area  compared  to  the  original  mass  estimate,  allows  information  regarding  site 
characterizations.  Site  age  considerations  and  plume  longevity  may  also  provide  insight  for  site 
assessments. 

DNAPL  sites  include  difficulties  and  uncertainties  pertaining  to  site  assessment.  In  order  to 
better  characterize  DNAPL  sites,  advancements  for  better,  faster,  and  cheaper  methods  for 
contaminant  spatial  distribution  and  integrated  spatial  or  temporal  measurements  have  been  the 
primary  focus.  The  first  method  offers  a  more  detailed  set  of  measurements,  but  becomes  costly. 
The  second  method  may  lack  the  spatial  or  temporal  resolution  necessary.  Cost  constraints  and 
the  necessity  of  quantity  of  data  determine  which  method  to  use,  resulting  in  a  tradeoff  from  the 
cost  effectiveness  and  the  amount  of  data  being  utilized. 
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The  quantity  and  quality  of  subsurface  site  data  compared  with  site  hydrogeology  and 
geochemistry  control  site  management  decisions  along  with  risk  assessment,  regulatory 
requirements,  and  long  term  goals. 

Prior  to  mass  flux  and  mass  discharge,  concentration  measurements  have  been  the  focus  of 
contaminant  regulations.  The  relationship  between  mass  flux  and  plume  response  is  the  key 
component  to  changing  the  focus  from  target  level  concentrations.  This  quantitative  relationship 
provides  a  better  understanding  of  the  site  and  allows  for  more  sound  decision  making.  This  also 
allows  for  a  better  understanding  of  the  flux  distribution  from  source  zones,  which  may  be  more 
cost  effective  for  remedial  designs  by  focusing  on  the  zones  that  provide  the  most  significant 
load  to  the  plume,  rather  than  the  areas  with  the  highest  concentration,  which  may  not  be 
contributing  significantly  to  the  plume.  Prioritizing  limited  resources  may  also  be  necessary  for 
sites  where  mass  cannot  be  removed  completely,  and  targeting  regions  responsible  for  plume 
behavior  must  be  the  focus. 


The  ITRC  DNAPL  team  identified  six  areas  where  mass  flux  and/or  mass  discharge  can 
improve  site  management  (ITRC,  2010): 

•  Site  conceptual  model  and  site  characterization 

•  Threat  and  exposure  evaluation 

•  Remediation  selection  and  design 

•  Performance  monitoring  and  optimization 

•  Compliance  and  long-term  monitoring 

•  Site  prioritization 

There  may  be  sites  with  multiple  improvements,  but  the  most  vital  is  the  understanding  of  the 
site. 

Understanding  DNAPL  sites  by  the  collection  of  contaminant  flux  data  allows  the  improvement 
of  site  management  and  decisions.  Flux  data  improves  the  ability  of  risk  assessment  and 
implementation  of  remedial  actions.  Traditional  site  data  collection  and  historical  site  data  can 
estimate  contaminant  mass  in  source  zones  and  plumes,  and  can  provide  a  better  understanding 
for  source  vs.  plume  remediation. 


b.  Mass  Discharge  and  Mass  Flux  as  Measures  in  Site 
Characterization  and  Management 

Mass  flux  and  mass  discharge  describe  the  transport  of  heat,  mass,  or  volumes  of  contaminant 
from  the  source  to  the  plume.  A  system  mass  balance  includes  those  terms,  and  defines  the  system 
conditions.  Mass  flux,  J  [M/L2/T],  is  defined  as  the  product  of  the  Darcy  flux,  g[L/T],  and  the  local 
concentration  of  contaminant,  C  [M/L3],  in  the  aqueous  phase 
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J=qC  . 


(Eq.  1) 


The  Darcy  flux  can  be  calculated  applying  Darcy’s  Law  along  the  direction  of  the  plume  axis 
at  the  location  of  interest: 


q  =  -K 


dh 

dx 


(Eq.  2) 


K  is  the  hydraulic  conductivity  [L/T],  and  h  is  the  head  measurement  at  the  point  of  interest.  The 
x  direction  is  selected  as  the  maximum  gradient  direction.  Darcy  flux  and  mass  flux  are  vector 
quantities,  and  thus  require  a  direction.  Typically  the  average  flux  vectors  are  associated  with  the 
flux  direction.  Although  groundwater  flow  is  dynamic,  and  contaminant  mass  flux  changes 
direction,  the  focus  is  on  a  large  frame  scale,  resulting  in  an  average  Darcy  flux  direction. 

The  contaminant  mass  flux  is  calculated  across  a  control  plane,  or  a  transect,  which  is 
perpendicular  to  the  plume  axis.  Integrating  the  mass  flux  with  the  transect  as  the  spatial  extent, 
quantifies  mass  discharge,  and  results  in  a  scalar  quantity.  A  positive  value  represents  the  transport 
in  the  plume  growth  direction. 


M 


D 


=  ^JdA 


(Eq.  3) 


Mass  discharge  may  vary  with  position  and  time.  Source  zone  mass  discharge  is  that  of  the 
down-gradient  edge  of  the  source  zone,  and  how  it  changes  with  natural  dissolution  results  in  the 
source  strength  function,  Mu(0,t).  The  relationship  between  the  source  zone  and  the  plume  can 
be  used  to  estimate  source/plume  mass  balance  which  in  turn  can  be  used  to  assess  site 
conditions. 

Site  age  is  again  defined  as  the  fraction  of  the  initial  source  zone  mass  that  has  been  removed 
from  the  source  zone  through  dissolution,  volatilization,  or  degradation,  and  may  contribute  to 
site  evaluation.  The  solubility  limit  of  the  contaminant,  the  groundwater  flow  velocity,  the  size  of 
the  source  zone  within  the  flow  direction,  and  the  average  NAPL  saturation  along  the  flow  path 
may  affect  the  site  age. 

The  average  of  the  mass  flux  over  sections  of  a  cross-section  may  be  useful  for  large  scale 
evaluations.  The  integral  pump  test  utilizes  multiple  wells  within  a  transect.  The  average  mass  flux 
over  an  area  B  can  represent  a  potential  impacted  well  from  a  capture  zone  of  the  same  area. 


\JdB 

7  =  ^ - 

B  (Eq.  4) 
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Average  concentrations  based  on  flux  can  be  calculated  as  C  =  J/q  which  relates  flux  to 
concentration,  which  may  be  useful  for  target  concentrations. 

Contaminant  distribution  among  phases  and  projection  of  mass  discharge  history  are  vital  for  site 
management.  Source  depletion  occurs  with  an  increasing  of  mass  discharge,  and  future 
projections  depend  on  mass  discharge  dynamics.  Declining  mass  discharge  suggests  that  the 
majority  of  the  mass  is  depleted  from  the  source  zone,  and  is  now  within  the  plume. 

Although  mass  discharge  rates  are  dynamic  and  attenuation  rates  also  vary,  hydrodynamic 
characteristics  of  the  aquifer  and  biogeochemical  conditions  also  contribute  to  the  dissolved 
plume  characteristics.  Steady-state  (no  growth)  conditions  occur  when  the  contaminant  mass 
discharge  from  the  source  zone  is  equal  to  the  overall  mass  attenuation  rate  within  the  plume.  If 
the  attenuation  is  greater  than  the  mass  discharge,  then  the  plume  is  expected  to  shrink,  but  if  the 
attenuation  rate  is  less  than  the  mass  discharge,  then  the  plume  is  expected  to  grow. 

Mass  discharge  history  may  be  obtained  from  measurements  of  the  source  zone  and  plume.  This 
history  provides  knowledge  about  the  site  and  offers  possible  projections  of  mass  discharge  of 
the  DNAPL. 

3.  Methods  for  Evaluating  DNAPL  Source  Dissolution 

a.  Relationship  between  Mass  Depletion,  Mass  Discharge  and  Mass  Flux 

Source  treatment  and  its  impact  on  mass  discharge  has  been  examined  with  modeling  analysis, 
specifically  the  REMCHLOR  model.  Due  to  hydrogeological  conditions  and  DNAPL 
distributions,  a  variety  of  changes  may  result.  After  removal  of  source  mass,  contaminant  mass 
discharge  reduces  significantly  according  isolated  testing  cells.  Actual  field  measurements  may 
have  alternate  results. 

i.  Power  Law  Model 

The  power  law  source  depletion  model  (PLM)  in  coordination  with  the  REMCHLOR  model, 
relates  the  flux-averaged  concentration  to  the  source  mass  empirically,  r  is  a  function  of  both 
media  and  DNAPL  variability. 


CfV) 

fcCs 


M(T ) 
Mn 


(Eq.  5) 


where  M(T)  is  the  mass  remaining  in  the  source  zone  [M]  and  Mo  is  the  initial  source  mass  [M]. 
Lalta  et  al  (2005a)  derived  an  analytical  solution  for  dissolution  under  the  PLM, 
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Cf(T) 

fcCs 


1 


(F  l)^D'0t  +  Ml-r 
Ml 


r 

l-r 


where  Md,o  is  the  initial  mass  discharge  (M/T). 


(Eq.  6) 


ii.  Streamtube  Model 

The  equilibrium  streamtube  model  (ESM)  is  an  analytical  model  proposed  to  simulate  DNAPL 
dissolution,  and  is  based  upon  a  Lagrangian  approach. 


fA 


1 . 1eJlnT-^) 


(Eq.  7) 


Here  C f  [ML"3]  is  the  flux-averaged  concentration  existing  the  source  zone;  T  is  the  number  of 
pore  volumes  of  natural  flowing  groundwater  that  have  moved  through  the  source  zone  [-];  fc  is 
the  fraction  of  the  streamtubes  initially  containing  DNAPL  [-],  Cs  is  the  solubility  limit  [ML'3]; 
t  is  defined  as  reactive  travel  time  that  combines  the  travel  time  (t)  and  the  trajectory-averaged 

$r/ 

DNAPL  content  5(5  =  ,v  ),  where  5V  is  the  trajectory-averaged  DNAPL  saturation  along  a 


streamline  [L3L3],  *>  is  the  porosity  [L3L3],  and  0W  is  the  water  content  [L3L  3];  //,„  r  and  crlri  r  arc 
the  mean  and  standard  deviation  of  the  log  normal  variable  In  r  .  The  travel  time  and  the  DNAPL 
content  5  are  both  assumed  to  be  lognormal  variables,  so  the  reactive  travel  time  obeys  the 
lognormal  distribution. 


iii.  Mass  Transfer  Model 

Parker  and  Park  (2004)  presented  a  simplified  model  for  estimating  DNAPL  source  zone 
mass  depletion  using  an  effective  Damkohler  number,  Da: 


CfV) 

cs 


l-exp[-Da] 


(Eq.  8) 


Here,  Da  =  keffLs  /  q  where  keff  [T1]  is  the  field-scale  effective  mass  transfer  coefficient,  Ls 

is  the  source  zone  length  in  the  mean  flow  direction  [L],  and  q  [LT"1]  is  the  average  Darcy  flux 
for  the  source  zone.  The  effective  mass  transfer  coefficient  ( keff )  is  related  to  changes  in  the 

DNAPL  mass  and  the  average  groundwater  velocity  (Parker  and  Park,  2004). 
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Figure  2.1.  Power  law  (PLM)  and  Equilibrium  Streamtube  (ESM)  models  for  source  zone 
dissolution. 


a.  Site  Mass  Balance  Considerations 

i.  Plume,  Source  and  Diffusion  Domain  Mass  Distribution 

Mass  from  the  plume  is  linked  with  mass  from  the  source  zone  through  the  historical  source 
strength  function.  Mass  flux  and  mass  discharge  relate  to  site  wide  mass  balance,  which  can  be 
useful  when  using  estimates  for  the  site. 

ii.  Concept  of  Site  Age  and  Aging 

Source  mass  estimation  can  be  estimated  using  an  exponential  source  decay  model  using  data 
from  wells  located  near  the  source. 


Mt2  =  Mn  exp(-(t2  -tx)k)  =  qA^,]  exp (~(t2  -  t{)k)  (Eq.  9) 

k 

where  Mti  and  Ma  are  the  source  masses  [M]  at  two  different  times,  and  Cti  is  the  concentration 
[ML-3]  in  the  monitoring  well  at  time  ti  [T].  Using  equation  7.9,  the  present  DNAPL  source  mass 
(Mti)  can  be  estimated. 


Another  exponential  source  decay  model  assumes  zero  degradation  rates  in  the  plume,  a  measure 
of  source  strength  at  some  time  ti ,  and  an  estimated  plume  mass  at  that  time. 


M 
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(Eq.  10) 


Md, 2  =MDn  eXP 
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(Eq.  11) 
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There  are  five  methods  discussed  that  approximate  source  strength  functions  at  field  sites,  and 
they  include  spatial  and  temporal  concentration  data. 

iii.  Mass  Discharge  and  Flux  during  Plume  Evolution 

Using  an  advection-dispersion  model  along  with  a  source  strength  function  provides  the  ability 
to  implement  site  characterization  from  field  data  that  was  collected  from  years  of  investigation. 

Wells  located  closer  to  the  source  zone  better  represent  source  zone  behavior  due  to  influence  on 
contaminant  concentration  in  areas  further  from  the  source  zone.  A  maximum  spatial  distance 
may  be  used  to  better  reflect  source  zone  behavior,  but  only  if  there  are  a  sufficient  number  of 
measurements. 

4.  Estimating  Source  Strength  Functions  at  Field  Sites 

Mass  discharge  can  vary  significantly  between  sites,  and  mass  flux  can  also  vary  within  the  same 
site.  Quantifying  source  strength  functions  is  discussed  in  the  next  sections. 

a.  Use  of  Time  Series  Data 

i.  Considering  Entire  Plume  Water  Quality  Data 

One  method  of  fitting  an  unknown  source  strength  function  to  the  concentration  time  series  field 
data  is  to  fit  the  data  globally,  and  try  to  optimize  the  parameters  r or  cr/nr.  With  this  method 
there  is  likely  a  distribution  of  parameters  that  must  be  used  to  represent  plume  response,  most 
often  at  higher  aged  sites  or  sites  with  a  large  contaminant  plume. 


ii.  Individual  Well  Water  Quality  Data  Series 

Contrary  to  global  optimization,  the  data  from  individual  wells  can  used  to  fit  source  strength 
function  parameters.  The  distribution  of  parameter  values  can  be  used  to  estimate  a  range  of 
values;  however  the  uncertainty  is  higher  due  to  outliers  from  other  processes  occurring  within 
the  plume.  Optimization  for  each  well  for  a  limited  number  of  parameters  can  be  more  effective 
and  less  computationally  intensive. 


b.  Use  of  Transects 

Wells  must  be  placed  perpendicular  to  the  flow  field  down-gradient  of  the  source  zone  when 
using  transects.  Transects  allow  mass  discharge  data  to  be  collected  temporally  or  spatially,  and 
contains  minimal  effects  of  noise  in  concentration  data.  Mass  discharge  across  transects  can  be 
calculated  using  a  number  of  methods,  but  in  general  it  can  be  summarized  as 

MD=fjq,C,Ai  .  (Eq.  12) 

1=1 
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Mass  fluxes  can  be  measured  using  passive  flux  meters,  multilevel  samplers,  stochastic  methods, 
or  a  combination  of  these.  Mass  discharge  can  be  calculated  by  interpolation  from  mass  flux  data. 
A  rigorous  stochastic  analysis  may  improve  accuracy. 

i.  Temporal  Transects 

Observing  the  change  in  flux  values  in  time  results  in  a  mass  discharge  time  series,  which  can  then 
be  fit  to  source  strength  functions.  One  assumption  being  made  is  that  the  transect  is  close  enough 
to  the  source  zone  to  reduce  error. 


ii.  Spatial  Transects 

Transect  data  collected  in  space  uses  multiple  well  transects  spaced  down  the  plume  centerline. 
This  allows  for  a  translation  of  concentrations  from  space  to  time,  and  is  done  using  travel  time 
using  assumed  retardation  and  groundwater  velocity  between  transects.  Losses  must  be  considered 
such  as  complete  degradation,  changes  in  hydro-stratigraphy,  and  remediation. 

c.  Use  of  Extraction  Wells 

Changes  in  concentrations  in  time  in  extracted  contaminated  groundwater  can  be  analyzed  for 
changes  in  source  strength,  which  can  predict  efficiencies  and  source  behavior  from  remediation. 
Although  it’s  possible  to  fit  a  source  strength  function  to  the  cumulative  mass  removal,  natural 
groundwater  gradients  may  not  be  reflected  by  the  source  strength  function  of  this  pumping. 

The  exponential  decay  model  is  a  simple  model  that  data  can  be  fit  to  represent  source  strength 
function  characteristics. 


M  A  _  J  _Cq 
M  a,o  Jo  Jo 


(Eq.  13) 


Here  Ma  is  the  mass  of  contaminant  per  unit  area  perpendicular  to  the  groundwater  Darcy  flux  q, 
J  is  the  mass  flux,  C  is  contaminant  concentration  and  the  0  subscripts  denote  initial  values.  Thus 
writing  an  equation  for  mass  loss: 

=  ~Cq  (Eq.  14) 

dt 


The  solution  for  the  exponential  decay  model  can  be  written: 


f 

C  -C0  exp 

v 


°—t 

Ma,0 


(Eq.  15) 


This  function  can  be  easily  fit  to  a  concentration  time  series  data  to  determine  the  unknown  values 
Co  and  (Jo/ Ma,o). 
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The  more  general  Power  source  function  model  can  also  be  applied  to  C(t)  data.  The 
solution  used  to  fit  to  C( t)  data  employs  the  following  relationship  between  source  zone  mass  and 
mass  flux: 


J_ 

X 


Ma 

Afi  J 


(Eq.  16) 


This  can  be  solved,  yielding 


c  =  c0 


J  o(r-i) 


r 

Er 


(Eq.  17) 


This  function  with  three  fitting  parameters  Co,  /"and  Jo/ Ma,o  can  also  be  fit  to  the  field  data 

sets. 


5.  Methods  for  mass  discharge  and  mass  flux  measurement 

Spatial  and  temporal  scales  of  measurement  are  two  methods  that  quantify  local  contaminant 
flux  and  integrated  mass  discharge.  Mass  flux  can  be  calculated  from  point  scale  techniques 
within  screened  wells,  and  that  data  can  be  spatially  integrated  to  determine  mass  discharge. 
Integral  pump  tests  collect  data  over  fully  screened  wells  and  large  volumes  of  water  pumped 
from  aquifers,  and  although  it  lacks  spatial  resolution,  “plume  cores”  are  incorporated. 


Point  Scale  Mass  Flux 


Integral  Scale  Mass  Flux 


Multilevel  Sampler  Passive  Flux  Integral  Pump  Test 

Method  Meter  Method  Method 


Figure  4.1  Point  and  integral  mass  flux  methods  in  wells  screens  typically  5  to  30 
feet  in  length  and  well  spacing  10  to  50  feet. 
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Local  measures  of  mass  flux  can  be  determined  from  local  flux  average  concentration,  local 
hydraulic  conductivity,  and  the  local  head  gradient  applied  in  the  maximum  flow  direction,  but 
typically  average  values  are  used  for  hydraulic  conductivity  and  the  head  gradient.  Spatial 
variability  is  obtained  from  the  concentration  distribution,  determined  using  multilevel  samplers. 


a.  Transect  Method 

The  transect  method  calculates  contaminant  discharge  using  the  product  of  the  well  cross  sectional 
area  within  a  transect  of  wells,  the  contaminant  concentration,  and  the  Darcy  flux.  The  Darcy  flux 
is  the  product  of  the  hydraulic  gradient  and  the  hydraulic  conductivity.  Average  values  are  used 
for  the  hydraulic  conductivity,  the  hydraulic  gradient,  and  the  flux  averaged  concentration  that  is 
collected  from  each  well.  The  contaminant  flux  is  the  product  of  the  Darcy  flux  and  the 
concentration  collected  from  each  well.  The  contaminant  flux  has  spatial  variability,  so  the  local 
value  may  be  significantly  different  from  the  average  value.  Integral  values  are  assumed  to  be 
acceptable  for  calculating  average  fluxes  and  integrated  fluxes  over  the  well  transect,  but  this 
assumption  has  not  been  tested  in  field  sites. 

The  transect  method  can  also  be  applied  to  areas  down  gradient  of  the  source  area  or  across  a 
plume,  using  the  local  mass  flux  and  local  contaminant  concentration.  The  Darcy  flux  can  only  be 
calculated  if  the  local  independent  values  are  measured.  The  local  hydraulic  conductivity  can  be 
measured  from  borehole  flow  meter  tests  or  slug  tests,  but  typically  average  values  are  used.  The 
local  hydraulic  gradient  is  rarely  determined,  so  average  values  are  used  from  monitoring  wells. 

The  well-established  methods  used  to  collect  data  offer  an  advantage.  The  transect  method  also 
provides  spatial  information,  and  it  incorporates  infrastructure  already  present  at  the  site  all  while 
producing  minimal  waste. 

The  transect  method’s  disadvantages  include  the  necessity  of  independent  measures  for  the  Darcy 
flux,  and  the  possible  errors  arising  from  sampling  small  volumes  of  groundwater,  and  in  order  to 
reduce  these  errors  more  measurements  must  be  taken.  Low  gradient  systems  can  introduce  errors 
for  independent  measures  of  flux  since  there  is  a  small  volume  of  groundwater  being  measured  at 
one  point  in  time,  so  data  must  be  integrated  to  obtain  contaminant  mass  loads.  A  study  by  Kubert 
estimated  a  range  of  20  to  50  percent  error  for  the  transect  method  when  the  “plume”  covers  a 
relatively  large  portion  of  the  control  plane.  A  recent  study  using  multilevel  transects  contained  a 
larger  margin  of  error  than  an  integrated  pump  test. 


b.  Passive  Flux  Meters 

Passive  Flux  Meters  (PFM)  use  a  sorbent  pack  that  covers  the  vertical  extent  of  the  contaminant 
plume,  and  they  allow  the  natural  gradient  groundwater  flow  during  a  designated  time  period  to 
capture  dissolved  contaminants.  At  the  time  of  placement,  PFM  are  pre-saturated  with  “resident 
tracers.”  During  the  process,  groundwater  fluxes  are  calculated  from  the  amount  of  depleted  tracer 
mass,  and  contaminant  fluxes  are  calculated  from  the  amount  of  contaminant  mass  captured.  The 
vertical  profile  of  the  PFM  allows  for  a  depth  distribution  to  be  determined. 
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Passive  Flux 
Meter 
Technology 


~  Retrieval  wire 
"  Tube  for  flow  bypass 

—  Sorbent  with  Tracers 
(activated  carbon) 

■  Viton  Washeis 
(minimize  vertical  flow) 


Figure  4.2.  Installation  of  the  passive  flux  meter  in  a  screened  well. 

Existing  site  conditions  allow  estimation  for  groundwater  flux  and  contaminant  concentrations, 
which  may  help  to  select  appropriate  tracers  so  that  tracer  loss  and  contaminant  mass  accumulation 
are  quantifiable. 

The  PFM  method  allows  spatial  information  to  be  collected,  integration  of  information  over  time 
can  be  done,  and  there  is  little  error  in  integrated  mass  load.  Based  on  measured  hydraulic 
gradients,  aquifer  conductivity  can  be  determined  from  this  method.  Minimal  waste  is  produced, 
and  this  method  can  be  applied  in  remote  settings. 

The  PFM  method  does  not  take  into  account  large  volumes  of  the  aquifer,  so  the  error  margin  may 
be  large.  The  data  must  be  integrated  in  order  to  determine  contaminant  mass  discharge  and  total 
water  discharge.  Tracers  that  are  required  must  be  approved  for  use  in  the  subsurface,  and  PFM 
have  not  been  applied  to  fractured  media  environments. 

This  approach  has  been  implicated  in  over  40  sites,  with  the  largest  at  Ft.  Fewis,  Washington, 
USA,  where  a  PFM  of  700  linear  feet  was  deployed  in  10-well  transects.  The  deepest  application 
took  place  near  Perth,  Australia  at  67  meters  deep. 


c.  Integral  Pump  Tests 

The  integral  pump  test  (IPT)  places  multiple  pumping  wells  perpendicular  to  the  direction  of 
groundwater  flows  and  measures  contaminant  concentration-time  series.  The  measured 
concentration  of  contaminants  is  used  with  the  natural  gradient  Darcy  flux  in  order  to  determine 
mass  flux  or  discharge.  The  significant  volume  of  water  allows  an  integral  measure  of  contaminant 
concentration  along  the  well  transect,  and  avoids  missing  “plume  cores”  within  the  transect. 

The  IPT  estimates  contaminant  mass  flow  as  well  as  supplies  information  on  the  spatial 
distribution  of  mass  flux  across  the  well  transect.  This  approach  can  be  applied  to  deep  aquifers, 
and  has  no  limitations  based  on  well  design. 
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The  IPT  generates  a  large  volume  of  waste  that  can  be  costly  depending  on  site  infrastructure. 
Tandem  circulation  wells  have  been  applied  in  some  sites  in  order  to  overcome  this  challenge.  The 
deployment  time  can  be  lengthy,  and  the  method  requires  independent  measures  of  the  Darcy  flux. 

6.  Ft.  Lewis  DNAPL  Case  Study  for  Mass  Flux  Evaluation 

The  data  collection  at  Ft.  Lewis  demonstrates  a  DNAPL  contaminated  site  with  extensive 
characterization  and  source  remediation  in  which  flux  assessments  were  used.  The  purpose  of  the 
data  is  to  look  beyond  remedial  performance  assessment  in  order  to  see  how  the  data  would  benefit 
remedial  design  and  site  management. 

The  Ft.  Lewis  TCE  contaminant  plume  is  3  km  long  with  a  source  zone  at  the  East  Gate  Disposal 
Yard  (EGDY).  Three  primary  contaminated  regions  were  identified  and  evaluated  for  a  variety  of 
remedial  technologies;  resistive  heating  was  selected.  Thermal  Remediation  Systems  (TRS)  was 
used  over  a  period  of  3  years,  starting  with  NAPL  Area  1,  the  most  upgradient  source.  The  focus 
of  remedial  assessment  using  flux  data  was  on  NAPL  Areas  1  and  3,  and  NAPL  Area  2  was  not 
evaluated  by  pre-  and  post-flux  measurements,  but  alternative  measurements  were  taken. 


Figure  4.3.  TCE  contaminant  plume  at  Ft.  Lewis  c.  2003. 
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Figure  4.4.  Ft.  Lewis  NAPL  treatment  areas. 

The  EGDY  source  areas  1  and  3  were  evaluated  in  order  to  quantify  mass  flux  and  integral  mass 
discharge.  A  single  transect  of  10  wells  were  located  immediately  down  gradient  of  treatment  zone 
1.  Since  zone  3  contained  mass  entering  from  up  gradient  sources,  a  transect  of  8  wells  was 
constructed  immediately  up  gradient  of  the  treatment  area.  Two  transects  of  10  wells  were 
constructed  immediately  down  gradient  of  zone  3  in  order  to  cover  the  saturated  thickness  of  the 
surficial  aquifer  of  50  ft. 

The  mass  discharge  was  evaluated  for  its’  effect  on  the  entire  plume  system.  A  site  wide  integral 
pump  test  was  conducted  immediately  down  gradient  of  the  source  area;  the  test  contained  15  fully 
screened  wells  and  was  conducted  just  after  conductive  heating  of  area  2.  The  measurements 
represent  an  intermediate  mass  discharge  value  that  was  impacted  by  remedial  efforts  in  Areas  1 
and  2. 

The  three  well  transects  contained  passive  flux  meters  followed  by  integral  pump  tests.  The  IPT 
initial  samples  represent  low  flow  samples  of  the  fully  screened  wells,  and  were  used  to  calculate 
mass  discharge  based  on  the  traditional  transect  method.  There  are  three  approaches  for 
quantifying  mass  discharge  for  each  flux  measurement,  and  each  approach  reports  values  for 
individual  wells  within  the  transect.  The  values  are  summed  to  quantify  mass  discharge. 


a.  Mass  Balance  Considerations 

Mass  distribution  within  the  source  area  and  within  the  plume  is  important  for  analyzing 
projections  of  mass  discharge  and  for  considering  short  term  and  long  term  site  management. 
Initial  mass  estimates  for  the  site  were  1 1 800,  13400,  43100  Kg  of  TCE  for  each  of  the  three  areas. 
Resistive  heating  removed  2580,  1090,  840  Kg,  which  demonstrated  that  the  initial  estimates  were 
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high  for  the  treatment  areas  due  to  the  fact  that  diminishing  returns  were  observed  for  the  mass 
recovery,  which  suggested  that  high  fractions  of  the  original  mass  were  removed.  Mass  removals 
during  resistive  heating  are  considered  to  be  more  reliable. 

Mass  discharge  from  NAPL  Areas  1  and  3  estimate  mass  discharge  history  for  each  of  the 
treatment  areas,  based  upon  certain  assumptions.  The  observed  mass  discharge  measured  at  NAPL 
Areas  1  and  3  were  6700/7300  mmoles/day,  nd  2300/8800  mmoles/day.  A  model  must  be  assumed, 
and  the  simplest  model  is  an  exponential  decay  model.  The  mass  at  the  beginning  of  the  treatment 
was  2600  KG,  and  the  mass  after  treatment  was  19  Kg,  resulting  in  a  removal  of  99.3%.  An  initial 
date  of  1950  was  assumed,  and  an  initial  mass  in  NAPL  Area  1  was  estimated  as  5200  Kg.  The 
same  analysis  for  Area  3  resulted  in  an  initial  mass  of  2800  Kg.  Based  upon  mass  removal,  NAPL 
Area  2  was  estimated  to  have  an  initial  source  mass  of  3100  Kg.  The  total  initial  estimate  of  10,000 
Kg  of  TCE  is  calculated,  and  approximately  5000  Kg  present  in  the  dissolved  plume. 

A  site  wide  IPT  conducted  6  months  after  heating  was  terminated  in  Area  1  provides  valuable 
insight  on  the  site-wide  mass  balance.  The  northern  and  central  sections  can  be  compared  with 
NAPL  Area  3.  The  site-wide  IPT  mass  discharge  of  3450  mmoles/day  compared  with  3200 
mmoles/day,  whereas  the  PFM  resulted  in  8800  mmoles/day.  NAPL  Area  1  must  be  summed  with 
an  unknown  estimate  from  NAPL  Area  2  to  estimate  the  site-wide  mass  discharge  prior  to 
treatment  since  NAPL  Area  2  is  down  gradient  from  NAPL  Area  1,  and  the  mass  discharge  may 
be  reduced  due  to  the  presence  of  TCE  in  solution  from  NAPL  Area  1.  Two  alternate  cases  may 
be  considered;  no  additional  mass  can  be  accounted  from  NAPL  Area  2,  or  NAPL  Area  2 
contributes  a  scaled  mass  input  based  on  results  from  other  zones.  These  two  cases  result  in  a  mass 
discharge  of  1 1000  to  14000  mmoles/day. 


b.  Remedial  Design  Considerations 

NAPL  Areas  1  and  3  were  evaluated  for  local  flux  measurements  using  transects  of  wells  using 
PFMs  segmented  in  vertical  intervals  of  30  cm.  At  the  down  gradient  section  of  NAPL  Area  1, 
each  well  contained  a  vertical  flux  profile  of  25  samples.  The  data  along  the  transect  was  used  to 
produce  a  contour  plot  of  the  flux  distribution,  and  the  data  demonstrated  the  variability  of  TCE 
mass  flux  emanating  from  the  DNAPL  source  area.  The  majority  of  the  mass  flux  is  localized 
within  the  control  plane,  but  there  are  likely  individual  pools  of  DNAPL  located  up  gradient.  In 
NAPL  Area  1,  about  90%  for  the  mass  flux  is  contained  within  about  33%  of  the  control  plane. 
This  demonstrated  the  mass  is  originating  up  gradient  of  the  flux  observation  at  the  control  plane, 
and  can  reduce  the  size  of  the  remediation  zone  by  focusing  on  the  up  gradient  regions. 
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Figure  4.5.  Ft.  Lewis  NAPL  treatment  area  1. 

The  control  plane  flux  data  is  considered  an  instantaneous  measurement  that  is  the  result  of 
groundwater  flow  imposed  on  the  DNAPL  source  area  for  the  travel  time  required  between  source 
and  control  plane,  taken  in  the  predominant  direction  of  groundwater  flow.  The  travel  time  for 
NAPL  Area  1  is  approximately  the  travel  time  through  the  treatment  area,  and  the  distance  is  about 
50  m.  PFMs  estimate  Darcy  velocity  of  0.25  m/day,  and  a  porosity  of  0.33,  resulting  in  a  travel 
time  of  70  days.  The  predominant  gradient  from  the  past  70  days  should  be  used  to  locate  up 
gradient  source  areas.  The  average  gradient  magnitude  and  direction  should  be  calculated,  but  this 
gradient  may  shift  between  wet  and  dry  seasons  so  this  must  be  considered  when  selecting  when 
to  collect  measurements. 

NAPL  Area  3  was  also  evaluated  both  up  gradient  and  down  gradient  for  flux  distributions,  and 
contour  maps  were  also  made.  The  flux  distribution  up  gradient  was  less  localized  than  Area  1 
most  likely  due  to  a  larger  travel  distance  between  the  source  and  the  control  plane  which  resulted 
in  more  mixing  and  dispersion. 
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Figure  4.6.  Shallow  down  gradient  mass  flux  distribution  NAPL  Area  3.  Flux  volumes  are 
mg/m2/d. 

The  down  gradient  control  plane  was  50  ft  in  saturated  thickness  which  was  much  thicker  than 
Area  1  which  is  the  result  of  the  till  unit  at  the  base  of  the  surficial  aquifer.  The  flux  distribution  is 
more  localized. 
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7.  Borden  Coal-Tar  Emplacement  Experiment  Case  Study  for  Mass 
Flux  Evaluation 

Chlorinated  solvent  DNAPL  sites  have  been  the  focus  of  mass  flux  evaluation;  however  the  same 
concepts  can  be  applied  to  manufactured  gas  plants  (MGP)  or  other  general  contaminated  sites 
with  a  source  and  a  plume  as  well  as  decisions  to  best  manage  the  site  with  limited  resources.  The 
relationship  between  the  source  strength  and  the  rate  of  degradation  represents  whether  a  plume  is 
advancing,  receding,  or  is  at  steady  state,  and  it  can  be  applied  to  sites  with  monitored  natural 
attenuation  of  polycyclic  aromatic  hydrocarbons  (PAH). 

There  are  differences  between  chlorinated  solvent  DNAPL  sites  and  coat  tar  DNAPL  sites, 
specifically  the  source  function  and  one  model  that’s  been  applied  is  the  power  law  model  (PLM). 
The  PLM  describes  source  function  behavior  with  a  screening -level  mathematical  description 
using  the  exponent,  T,  which  is  a  function  of  the  DNAPL  spatial  distribution,  the  hydraulic 
conductivity  spatial  distribution,  and  the  correlation  between  the  two.  Dissolution  behavior 
associated  with  multi-component  NAPL  is  not  accounted  for  in  this  model. 

Results  from  a  field  experiment  in  August  1991  where  74  Kg  of  coal  tar  creosote  was  injected  in 
a  shallow  aquifer  at  the  Canadian  Forces  Base  in  Borden,  Canada  will  be  discussed.  The  source, 
plume,  and  recent  results  have  been  monitored.  The  emplacement  involved  mixing  creosote  above 
grade  with  -5800  kg  of  coarse  sand,  then  the  material  was  placed  below  the  water  table  in  two  2m 
by  2m  by  0.5  m  blocks  separated  by  lm.  The  bulk  saturation  was  about  7%. 

It  was  assumed  the  DNAPL  was  uniformly  distributed  and  had  homogeneous  hydraulic 
conductivity,  which  would  theoretically  create  a  source  function  in  which  the  majority  of  the 
DNAPL  would  have  to  be  depleted  before  reductions  in  source  mass  discharge  occurred.  This 
behavior  represents  a  small  PLM  exponent.  Individual  components  in  a  multi-component  DNAPL 
behave  very  differently  even  though  a  multi-component  DNAPL  as  a  whole  behaves  similar  to  a 
single  component  DNAPL. 

One  relationship  used  to  describe  dissolution  from  a  multi-component  DNAPL  is  Raoult’s  law 
(Banerjee,  1984),  which  states  that  dissolution  of  the  ith  species  is  a  function  of  the  aqueous 
solubility  and  DNAPL  mole  fraction: 

C,  =  Z./.S,  ,  (Eq.  18) 

where  G  is  the  aqueous  concentration  [ML"3]  of  the  ith  solute  resulting  from  equilibrium 
partitioning  from  the  DNAPL;  and  yi  are  the  mole  fraction  [-]  and  the  activity  coefficient  [-] 
of  the  ith  solute  in  the  organic  phase,  respectively;  and  Si  is  the  pure-phase  aqueous  solubility  limit 
[ML'3]  of  the  ith  solute.  Often  yt  is  assumed  to  be  1,  which  is  the  value  in  an  ideal  organic  fluid. 

Raoult’s  law  has  been  applied  to  predict  solute  partitioning  between  coal  tar  and  water. 
Measurements  and  estimates  based  on  Raoult’s  law  have  been  in  agreement,  with  the  difference 
ranging  from  factors  of  2  to  4;  others  have  noticed  order-of  magnitude  differences  or  have  noticed 
kinetic  limitations  in  mass  transfer  between  water  and  coal  tar.  Another  method  to  predict  solute 
partitioning  between  water  and  coal  tar  is  polyparameter  linear  free  energy  relationships  which 
uses  a  mechanistic  description  of  the  intermolecular  interactions. 
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Raoult’s  law  has  been  applied  to  the  Borden  emplacement  study.  The  PLM  can  be  applied  for  the 
source  function  using  an  exponent  of  1  for  more  soluble  compounds,  but  the  less  soluble 
compounds  in  the  early  time  data  deviate  from  the  PLM  behavior  due  to  the  mole  fraction 
increasing  with  time  as  more  soluble  compounds  leave  the  DNAPL.  Any  less  soluble  compound 
becomes  more  soluble  after  some  time,  and  its  dissolution  follows  a  PLM  with  an  exponent  of 
unity. 
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Figure  6.1.  a)  Source  strength  function  (relative  mass  discharge  versus  time)  for  select  PAH 
compounds  based  on  Raoult’s  law  and  data  from  the  Borden  coal  tar  emplacement 
experiment,  and  b)  source  function  (relative  mass  discharge  versus  relative  mass)  for  the  same 
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select  PAH  compounds  based  on  Raoult’s  law  and  data  from  the  Borden  coal  tar  emplacement 
experiment  (King  and  Barker,  1999). 

Shown  in  Figure  7.10  is  mass  discharge,  Md,  as  a  function  of  time  for  10  compounds  as 
reported  by  Fraser  et  al  (2008);  6  measurements  were  completed  between  278  to  514  days 
after  emplacement.  This  data  provides  a  useful  opportunity  to  test  how  well  simple  models 
of  source  strength  functions  predict  future  behavior  using  intermediate  scale  pilot  tests.  It 
will  be  assumed  that  the  first  three  measurements  collected  between  278  and  1357  days  after 
emplacement,  were  part  of  site  monitoring  activities  and  will  be  used  to  make  predictions  of 
mass  discharge.  These  predictions  in  turn  are  compared  to  measurements  completed  for 
3619<f  <5140  days. 


Figure  6.2.  Source  mass  discharge  (g/day)  as  a  function  of  time  (days)  for  10  monitored 
compounds.  Data  reported  by  Frasier  et  al.  (2008). 


If  the  first  three  measurements  are  viewed  in  isolation,  it  would  be  assessed  that 
monotonically  decreasing  trends  in  Md  are  evident  for  m- xylene  and  naphthalene,  while 
monotonically  increasing  trends  in  Md  are  evident  for  phenanthrene  and  anthracene.  With 
the  exception  of  fluorine,  mass  discharge  for  the  remaining  compounds  increase  over  the  first 
two  measurements,  but  decrease  with  the  third.  Results  for  fluorine  were  not  reported  for  t  = 
1357  days,  but  Md  increased  over  the  first  two  measurements.  This  behavior  can  be  explained 
by  the  more  soluble  compounds  leaving  the  NAPL  which  causes  the  mole  ratio  of  the  less 
soluble  compounds  to  increase,  resulting  in  a  higher  effective  solubility.  It  is  noted  that  the 
response  at  the  source  zone  control  plane  is  an  integrated  measurement  process  that  represents 
dissolution  as  well  as  sorption/desorption  processes  between  NAPL  and  the  control  plane, 
and  potential  mass  loss. 
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M-xylene  and  naphthalene  show  decreasing  trends  that  can  be  fit  to  the  power  law  model. 
Model  fits  were  completed  for  F  =  {0.5,  1.0,  2.0}  using  linear  regression,  and  are  shown  in 
Figure  7.11.  For  F  =  1.0,  a  linear  transformation  of  the  mass  discharge  data  is  required  prior 
to  linear  regression  using  the  natural  logarithm  (i.e.,  y  =  In [Md  (t)],  where  y  is  the  transformed 
data  upon  which  linear  regression  is  completed).  Likewise,  for  T  =  2.0,  a  linear 
transformation  of  the  mass  discharge  data  is  needed  using  the  inverse  of  the  square  root  (i.e., 
y  =  Mo(t)"1/2),  where  y  is  the  transformed  data  upon  which  linear  regression  is  completed). 
As  shown  in  Figure  6.3,  the  mass  discharge  measurements  competed  for  3619  <t<  5140 
days  fall  in  between  model  fits  for  T  =  { 1.0,  2.0}.  The  model  fit  using  T  =  0.5  results  in  the 
largest  discrepancy  between  measured  data  and  model  predictions.  However,  the  rapid 
decline  in  the  source  strength  as  predicted  by  a  PLM  with  F  =  0.5  is  generally  considered 
unlikely. 


Figure  6.3.  PLM  fits  to  measured  data  for  m-xylene  and  naphthalene. 
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a.  Source  Mass  Estimation 

Integration  of  core  data,  partitioning  tracer  tests,  and  geophysical  methods  have  been  used  to 
estimate  source  mass;  however  estimates  have  high  uncertainty  due  to  spatial  heterogeneity  in 
subsurface  characteristics  and  contaminants. 

Once  a  PLM  has  been  fit  to  the  data,  it  can  be  used  to  estimate  the  initial  source  mass,  Mo  (i.e. 
Eq.  7).  For  m- xylene,  estimates  of  Mo  were  {2170,  2204,  3441 }  g  for  T  ={0.5,  1,2},  respectively; 
while  the  emplaced  mass  of  m- xylene  was  2850  g.  Estimates  of  the  initial  mass  were  therefore 
within  a  factor  of  2  of  the  true  mass.  For  naphthalene,  estimates  of  Mo  were  {3298,  4176,  6076} 
g  for  r  ={0.5,  1,  2},  respectively;  while  the  emplaced  mass  of  naphthalene  was  71 10  g.  Estimates 
of  initial  mass  were  therefore  within  a  factor  of  2.2  of  the  true  mass.  For  any  field-based  method 
of  DNAPL  mass  estimate,  these  results  have  to  be  considered  as  very  good.  However,  the 
controlled  nature  of  the  study  may  be  responsible  for  the  agreement  between  the  estimates  and  the 
actual  values,  and  applications  of  this  approach  to  more  complex  field  sites  need  further  evaluation. 

b.  Borden  Coal-Tar  Plume  Behavior 

Evaluation  of  the  plume  behavior  with  respect  to  source  strength  and  the  attenuation  capacity 
within  the  plume  was  important  for  flux  based  site  management.  Plume  behavior  for  the  first  four 
years  of  the  study  was  reported  by  King  and  Barker  (1999)  and  King  et  al.  (1999)  for  seven 
compounds:  phenol,  m-xylene,  naphthalene,  phenanthrene,  dibenzonfuran,  carbazole,  and  1- 
methylnaphthalene.  Phenol  (reported  solubility  of  82,000  mg/L)  was  observed  to  completely 
dissolve  from  the  source  and  migrate  downgradient  as  a  discrete  slug,  and  was  completely 
degraded  after  2  years.  For  m-xylene  and  carbazole,  the  plume  was  observed  to  migrate 
downgradient  and  then  recede  as  the  mass  flux  from  the  source  decreased  below  the  natural 
attenuation  capacity  within  the  plume.  The  spatial  extent  of  dibenzofuran  reached  steady  state, 
indicating  a  balance  between  the  mass  discharge  from  the  source  and  the  rate  of  plume  degradation. 
Napthalene  and  1-methylnaphthalene  continued  to  migrate  downgradient  over  the  entire  study 
period,  indicating  source  mass  discharge  exceeded  degradation.  The  phenanthrene  plume  was  also 
observed  to  expand,  and  recede,  however  the  degree  of  sorption  associated  with  this  compound 
hampered  assessment  over  the  study  duration.  Fraser  et  al.  (2008)  extended  the  analysis  presented 
by  King  and  Barker  (1999)  and  King  et  al.  (1999)  to  a  duration  of  14  years,  and  included  additional 
compounds  in  their  analysis  (anthracene,  acenaphthene,  biphenyl,  and  fluorine).  The  compounds 
were  placed  into  one  of  three  classes  based  on  their  plume  behavior.  The  first  class  consisted  of 
those  compounds  that  were  significantly  attenuated  over  the  14-year  study  duration  (phenol,  m- 
xylene,  and  naphthalene),  and  represents  the  expected  result  when  the  natural  attenuation  capacity 
exceeds  the  source  strength.  The  second  class  included  the  compounds  for  which  the  plume 
continued  to  expand  over  the  entire  14-year  study  (phenanthrene  and  dibenzofuran),  and  represents 
the  expected  result  when  the  source  strength  exceeds  the  downgradient  natural  attenuation 
capacity.  The  final  class  were  those  compounds  that  appeared  to  reach  a  steady  plume  extent  (1- 
methylnaphthalene,  carbazole,  biphenyl,  fluorine,  anthracene,  and  acenaphthene),  and  represents 
the  expected  result  when  the  natural  attenuation  capacity  and  source  strength  are  balanced. 

From  the  perspective  of  flux-based  site  management,  and  its  application  to  future  MPG  sites, 
the  Borden  emplacement  study  highlights  the  dependence  of  the  plume  and  its  spatial  extent  on 
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the  source  mass,  source  age,  source  strength,  and  natural  attenuation  capacity.  Manipulation  of 
source  mass  and  strength  (through  source  remediation),  or  plume  degradation  (through  plume 
remediation)  are  the  options  that  can  be  used  to  optimize  the  management  of  the  contaminated  site 
if  complete  DNAPL  source  depletion  is  not  feasible. 


8.  Alameda  Site  Case  Study  Using  Source  zone  Capture  Wells  to 
Enhance  DNAPL  Dissolution 


a.  Short  Term  Pumping  to  Estimate  Source  Strength  Function 

The  Alameda  Site  on  Alameda  Island  located  near  San  Francisco,  California  is  a  small,  shallow 
DNAPL  site  contaminated  with  VOCs,  mostly  TCE,  in  the  shallow  aquifer.  Research  studies 
were  conducted  to  create  a  conceptualized  site  model  and  to  determine  the  source  strength 
function  (SSF)  based  on  contaminant  dissolution  (see  Wang  et  al.,  2014).  The  study  used 
historical  groundwater  data  and  membrane  interface  probes  to  provide  an  initial  estimate  of  the 
extent  of  the  DNAPL  source.  Soil  cores,  photoionization  detectors  (PID),  and  a  hydraulic 
profiling  tool  (HPT)  were  also  used  to  determine  site  characterization.  The  data  showed  the  site 
contained  a  low  flow  zone  from  4.6  m  to  6.1  m,  and  a  high  flow  zone  from  6.1  m  to  7.6  m. 

Passive  flux  meters  (PFMs)  were  installed  into  wills  with  screened  intervals  in  both  the  low  flow 
and  high  flow  zone  in  the  source  zone  in  June  2010  with  a  vertical  sampling  resolution  of  15  cm. 
MLS  wells  were  also  installed  with  a  vertical  sampling  interval  of  0.30  m,  focusing  on  the  low 
flow  zone.  A  push-pull  tracer  test  was  conduction  in  August  2010  in  two  wells  within  the  source 
zone,  and  one  well  in  the  plume. 


Temporal  DNAPL  dissolution,  also  known  as  SSF,  can  help  determine  mass  flux  and  mass 
discharge  which  can  be  applied  to  long  term  effects  of  source  depletion  efforts.  The  equilibrium 
streamtube  model  (EST)  and  the  power  law  model  (PLM)  were  two  models  that  were  applied  to 
this  site.  The  EST  model  is  based  on  a  Lagrangian  approach  where  the  DNAPL  source  zone  is 
conceptualized  as  a  collection  of  non-interacting  streamtubes,  with  hydrodynamic  and  DNAPL 
heterogeneity  represented  by  the  variation  of  the  travel  time  and  DNAPL  saturation  among  the 
streamtubes  (Jawitz  et  al.,  2003;  Jawitz  et  al.,  2005). 


c(t)=\f,c, 


1-erf 


z,  \ 

ln^-Amr 


cr,, 


(Eq.  19) 


where  C(t )  is  the  flux- averaged  concentration  leaving  the  source  zone;  fc  is  the  fraction  of  the 


streamtubes  initially  containing  DNAPL;  Cs  is  the  solubility  limit  of  the  DNAPL  in  water;  T  is 
defined  as  reactive  travel  time  that  combines  the  travel  time  t  and  the  trajectory-averaged  DNAPL 


content  £  (S 


,  where  SN  is  the  trajectory-averaged  DNAPL  saturation  along  a  streamline, 


V  is  the  porosity  ,  and  9W  is  the  water  content);  //lnr  and  crlnr  are  the  mean  and  standard  deviation 
of  the  log  normal  variable  In  r  .  The  travel  time  t  and  the  DNAPL  content  S  are  both  assumed  to 
be  lognormal  variables,  so  the  reactive  travel  time  also  follows  the  lognormal  distribution. 
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The  EST  model  is  an  accurate  method  of  predicting  DNAPL  dissolution  by  using  partitioning 
tracer  tests  (PTT)  to  independently  parameterize  the  model.  The  initial  aqueous  concentration  of 
DNAPL  that  is  used  to  estimate  fc  by  dividing  by  the  solubility  limit  and  the  nonpartitioning 

and  partitioning  tracers  that  show  adequate  separation  are  the  data  needed  to  parameterize 
Equation  1.  Other  tests  such  as  soil  cores,  multilevel  samplers  (MLSs),  passive  flux  meters 
(PFM)  and  push-pull  tracer  tests  may  be  used  to  supplement  the  tracer  data. 


The  PLM  has  been  presented  by  Falta  et  al.  (2005b),  and  Zhu  and  Sykes  (2004), 

r 

|Er 


c(t)=f<c, 


(r  i )f,c,vaA 


M. 


t  +  1 


(Eq.  20) 


where  T  is  a  PLM  fitting  parameter;  Vd  is  the  Darcy  velocity;  A  is  the  cross-sectional  area  through 


which  contaminants  migrate;  and  M0  is  initial  mass  in  the  source  zone.  The  special  case  of  V  =  I 
leads  to  a  simple  exponential  form  for  the  source  zone  concentration  over  time  (Zhu  and  Sykes, 
2004), 


C(f)  =  /cC,exp 


fcCsVdA1 


(Eq.  21) 


In  laboratory  experiments,  Chen  and  Jawitz  (2009)  found  that  dissolution  follows  an  exponential 
decay  at  aged  contaminated  sites  where  a  significant  fraction  of  the  initial  DNAPL  mass  has  been 
eluted  from  the  source  zone. 


The  volume  of  contaminant  is  difficult  to  know  at  DNAPL  contaminated  sites;  however  an 
estimate  is  essential  to  accurately  estimate  SSF  parameters.  Source  mass  was  estimated  using  four 
methods  (M1-M4).  The  first  method  used  soil  core  data  combined  with  membrane  interface  probe 
data  (MIP)  to  delineate  the  extent  of  the  DNAPL.  The  soil  core  data  resulted  in  a  high  concentration 
of  TCE  (2200  mg/kg)  at  5.8  m  below  ground  surface  (bgs),  and  (6600  mg/kg)  at  6.1  m  bgs.  All 
other  core  samples  were  less  than  100  mg/kg,  which  resulted  in  no  DNAPL  present.  An 
approximate  minimum  (yellow  circle  area  in  Figure  7.1)  and  maximum  (square  area  in  Figure  7.1) 
extent  of  the  DNAPL  mass  was  delineated  and  the  mass  volume  was  determined  using  Equation 
22, 

VN  =  PbVsoilC  /  pTCE  (Eq.  22) 


where  ph  is  soil  bulk  density  and  pTCE  is  TCE  density  ;  Vsoil  is  volume  of  soil;  C  is  average 
concentration  based  on  the  high  resolution  soil  core  at  SPW-3-2  (mg/kg). 


The  second  method  used  partitioning  tracer  test  (PTT)  data  to  estimate  the  average  TCE 
concentration  in  the  water  swept  volume  based  on  the  separation  between  two  partitioning  tracers 
that  had  different  partitioning  coefficients.  The  average  saturation  of  TCE  in  the  source  zone  can 
be  estimated  by  the  first  normalized  moment  of  the  two  partitioning  tracers, 


SN  - 


R- 1 


(Eq.  23) 


where  the  retardation  factor,  R  ,  is  the  ratio  of  the  first  normalized  moment  of  the  partitioning 
tracers,  KN  is  the  tracer  TCE-water  partitioning  coefficient.  The  total  DNAPL  volume,  VN  , 
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within  the  swept  volume  of  each  well  was  determined  by  multiplying  the  average  SN  and  swept 
volume  Vs ,  determined  for  each  well  based  on  the  first  moment  of  n-hexanol  tracer  data  which 

has  a  lower  KN  to  avoid  the  high  uncertainties  and  complications  using  the  recirculated 
methanol  or  IP  A. 


The  third  method  used  both  tracer  test  and  aqueous  dissolution  data  during  the  push-pull 
tracer  test  and  use  of  EST  model  theory  to  estimate  DNAPL  mass.  The  separation  between  the 
partitioning  and  nonpartitioning  tracers  for  the  push-pull  tracer  test  was  not  distinct  enough  to 
directly  estimate  the  DNAPL  mass.  An  alternative  approach  was  applied  where  the  push-pull 
nonpartitioning  tracer  data  (IPA)  and  contaminant  dissolution  data  were  scaled  in  order  to  make 
the  solutes  have  similar  radial  travel  distance  and  comparable  relative  concentration  (Equation  6 
and  7). 

Cnp=l-^~  (Eq.  24) 

^0 -np 

C=-^~  (Eq.  25) 

peak 

where  Cnp  is  the  normalized  nonpartitioning  tracer  concentration;  C  is  nonpartitioning  tracer 
data  during  the  push-pull  tracer  test;  C0_np  is  input  tracer  concentration;  C  is  TCE  concentration 
C  scaled  to  the  peak  concentration  C  k  during  the  push-pull  tracer  test.  The  separation 

between  scaled  nonpartitioning  tracer  data  and  dissolution  data  indicates  the  presence  of  DNAPL 
in  the  tracer  swept  zone  and  potentially  can  be  used  to  estimate  the  mass  in  the  source  zone.  The 
procedure  proceeded  by  first  fitting  the  scaled  nonpartitioning  tracer  data  to  Equation  1  to 
determine  juUu  and  <rln,  by  assuming  no  DNAPL  (juXnr  =  //ln; ,  cr|nr  =  crlnr),  and  then  calculating 


the  first  moment  of  the  t  distribution  m[  using, 


ml  -  exp 


Mint  +- 


r2  A 

;  Inr 


(Eq.  26) 


The  peak  scaled  TCE  dissolution  data  was  fit  to  Equation  1  to  determine  /jlnT ,  alnr  ,  fc  and  then 


calculated  the  first  moment  of  the  r  distribution  m{  using. 


ml  =exp 


Mint  +' 


r2  A 

'  lnr 


(Eq.  27) 


A  s 

Then,  Equation  9  was  used  to  calculated  the  first  moment  of  the  S  distribution  m,  (Basu  et  al., 
2008a;  Jawitz  et  al.,  2003), 

m,r  =  Kwm[ml  (Eq.  28) 

where  Kw  is  equal  pTCE  /  Cs .  The  swept  volume  of  the  push-pull  test  was  smaller  than  the  source 

zone  volume,  but  in  this  estimation,  the  same  TCE  distribution  was  assumed  outside  and  within 
the  push  pull  swept  volume  as  in  the  source  zone.  Finally,  the  volume  of  DNAPL  in  the  source 
zone  was  estimated  using, 

VN  =  fcVSz>ni  (Eq.  29) 
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where  Vsz  is  the  source  zone  volume  (square  in  Figure  7.1)  which  is  different  from  the  swept 
volume  Vs . 

The  fourth  method  used  tracer  peak  concentration  data  during  the  PTT  collected  at  the 
MLSs  and  also  applied  to  the  extraction  wells.  Here,  a  simple  linear  partitioning  isotherm  was 
assumed  to  partition  tracer  between  the  volume  of  water  and  volume  of  DNAPL  in  the  MLS  flow 
path  as  if  the  system  was  a  batch  reactor.  Based  on  the  partitioning  relationship  between  NAPL 
and  water  phase,  and  mass  balance  of  nonpartioning  tracer  (methanol)  and  partitioning  tracer  (2- 
octonal),  NAPL  volume  was  derived, 

K„=^-  (Eq.30) 

^  p-w 

M  p-T  =  C0-pVin  =  VwCp-w  +  VNCp-N 

M  np-T  =  Co-llpV;n  =  Cnp_wVw 

Combining  these  equations  and  rearranging  produces 

C  C 

np-w  p-w 

v„  =  C°-’f  C°^  vw  =  S„VW  (Eq.  31) 

'•'p-w 

^NW  ^ 

^0 -p 

where  C p_N  and  Cp_w  are  tracer  concentrations  in  the  DNAPL  and  aqueous  phase;  M  p_T  is  the 
total  partitioning  tracer  mass;  M  T  is  total  nonpartitioning  tracer  mass;  Vw  and  Vjn  are  volume 
of  water,  volume  of  DNAPL  and  tracer  input  volume,  respectively. 

Ml  estimated  volume  at  a  range  of  2  L  -  20  L,  which  is  a  wide  but  reasonable  range.  M2  estimated 
11  L  overall,  but  had  a  retardation  factor  of  1.06  which  is  low  because  the  overall  saturation  was 
low.  The  swept  volume  was  49  m3.  The  highest  volume  estimated  was  5.8  L  with  a  retardation 
factor  of  1.12,  but  contains  the  highest  uncertainties  due  to  the  retardation  factor  being  below  1.2. 
M4-2  estimated  37  L  of  DNAPL,  which  is  much  greater  than  the  9.9  L  estimated  using  the  same 
method  applied  to  the  tracer  data  in  the  extraction  wells  (M4-1).  M4-1  could  have  underestimated, 
or  M4-2  overestimated  the  localized  mass  which  produced  larger  peak  differences.  This  method 
may  not  be  representative  of  the  entire  source  zone,  and  without  the  MLS  PTT  data,  the  volume 
of  DNAPL  in  the  low  flow  zone  may  not  have  been  detected  by  other  integral  methods. 

The  site  was  estimated  to  be  aged  greater  than  80%  based  on  site  mass  balance.  The  initial  mass 
was  estimated  using  the  methods  presented,  and  using  T  =  1  as  a  good  approximation.  The  EST 
used  the  parameters  estimated  by  the  PTT,  and  by  fitting  parameters  using  the  equations  presented 
and  the  push-pull  tracer  test.  The  aquifer  and  the  TCE  heterogeneity  in  the  swept  volume  was 
assumed  to  represent  the  entire  source  area. 
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b.  Use  of  Flux  and  Tracer  Data  to  Improve  Mass  Balance 

Predicting  aqueous  dissolution  requires  accurate  mass  in  the  source  zone  and  in  the  fraction 
contaminated  zone.  The  volume  of  contaminant  is  also  required  to  predict  SSF  parameters.  One 
method  to  calculate  DNAPL  volume  is  by  using  partitioning  tracers.  Particularly,  the  equilibrium 
streamtube  model  (EST)  uses  tracer  tests  to  estimate  parameters.  Partitioning  tracer  tests,  passive 
flux  meters,  multilevel  samplers  (MLS),  and  push-pull  tracer  tests  may  all  be  useful.  The  data 
collected  from  these  tests  during  the  early  phase  of  site  investigation  provide  essential  data  that 
can  help  delineate  and  confirm  the  extent  of  the  source  zone,  and  it  can  also  be  used  to  estimate 
source  mass.  All  of  the  data  combined  can  improve  parameterization  for  the  EST  model. 

Nonpartitioning  tracer  data  can  be  used  to  estimate  average  travel  time,  travel  time  distribution, 
and  to  calculate  the  swept  volume.  Average  pore  velocities  were  estimated  based  on  the 
partitioning  tracer  test  nonpartitioning  tracer  (methanol)  peak  arrival  time  at  MLSs.  As  discussed, 
the  data  demonstrated  low  flow  zones,  and  high  flow  zones  separated  at  a  depth  of  5.9  m  below 
ground  surface.  Passive  flux  meters  were  installed  into  wells  SPW-1,  SPW-2,  SPW-3,  and  SPW- 
4  with  screen  intervals  in  both  the  low-flow  zone  and  the  high-flow  zone.  Multilevel  sampler  wells 
were  installed  at  locations  SMLS-1  to  SMLS-8  with  a  vertical  sampling  interval  0.30  m,  focused 
on  the  low-flow  zone.  Passive  flux  meters  were  also  deployed  in  the  four  extraction  wells  six  and 
a  half  months  after  the  end  of  the  enhanced  dissolution  test.  The  velocity  of  the  contaminant  flow, 
as  well  as  the  source  mass,  allow  more  information  to  be  applied  to  mass  balance  of  the  site. 

Table  1.  Well  screen  intervals  and  radius  information  in  the  source  zone 


Well 

up 

Screen  Interval  (m  bgs) 

bottom 

Diameter 

(cm) 

SIM  1~3 

4.6 

7.6 

3.8 

SMLS  1~8 

5.0 

7.0 

1.3 

PEW  1-4 

4.6 

7.6 

10 

SPW-1  ~2 

5.2 

6.7 

3.8 

SPW-3- 1 

5.2 

5.9 

3.8 

SPW-3-2 

5.9 

6.7 

3.8 

SPW-4 

5.2 

6.7 

3.8 
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Table  2.  Tracer  partitioning  coefficient  Kn  (TCE  -water)  and  input  concentration  Co  for  the  push-pull 


tracer  test  and  partitioning  tracer  test  (PTT). 


Tracers 

Knw 

Co(mg/L) 

methanol 

0 

530 

IPA 

0 

1050 

Push  pull 

n-hexanol 

21 

470 

tracer  test 

DMP 

63 

540 

e-HEX 

280 

130 

1-octonal 

421 

130 

methanol 

0 

430 

Partitioning 

IPA 

0 

1000 

tracer  test 

n-hexanol 

99 

380 

(PTT) 

2-octonal 

266 

120 

e-HEX 

280 

120 
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Table  3.  NAPL  volume  estimates  using  four  methods  including  swept  volume,  NAPL  saturation  and  total 
mass  volumes. 


Vsoii  (m3) 

Vn  (L) 

Ml 

min 

4.0 

2.0 

max 

40 

20 

Vis  (m3) 

R 

sN 

Vn  (L) 

PEW-01 

15 

1.05 

1.99E-04 

3.0 

M2 

PEW-02 

12 

1.12 

4.74E-04 

5.8 

PEW-03 

10 

1.06 

2.48E-04 

2.5 

PEW-04 

12 

1.00 

5.27E-06 

0.06 

PEW-T 

49 

1.06 

2.32E-04 

11 

M3 

fc 

s 

mx 

Vs  (m3) 

Vn  (L) 

SPW-3-1 

0.14 

2.11E-03 

12 

3.7 

Vs  (m3) 

SN 

Vjv  (L) 

PEW-01 

15 

1.38E-04 

2.1 

M4-1 

PEW-02 

12 

4.49E-04 

5.4 

PEW-03 

10 

2.06E-04 

2.1 

PEW-04 

12 

2.92E-05 

0.35 

sum 

49 

9.9 

SMLS-4-1 

1.15 

0 

0 

SMLS-4-2 

1.15 

4.25E-03 

4.9 

SMLS-4-3 

1.15 

9.35E-03 

11 

SMLS-4-4 

1.15 

2.86E-04 

0.33 

SMLS-7-1 

0.60 

7.33E-03 

4.4 

SMLS-7-2 

0.60 

9.24E-03 

5.5 

M4-2 

SMLS-7-3 

0.60 

6.08E-03 

3.6 

SMLS-7-4 

0.60 

0 

0 

SMLS-8-1 

0.30 

9.96E-03 

3.0 

SMLS-8-2 

0.30 

9.37E-03 

2.8 

SMLS-8-3 

0.30 

6.81E-03 

2.0 

SMLS-8-4 

0.30 

0 

0 

sum 

8.17 

37 
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Table  4.  EST  model  parameters  and  volume  of  mass  in  the  swept  volume  for  each  individual  extraction 
well  and  flux  averaged  for  all  wells.  The  results  are  based  on  the  best  fit  of  the  enhanced 
dissolution  data  (mass  discharge)  with  the  EST  model  using  a  generic  algorithm  code. 


E 

fc 

&lnz 

Vn(L) 

PEW-l 

0.93 

0.016 

3.50 

1.32 

2.6 

PEW-2 

0.80 

0.024 

3.85 

1.33 

5.5 

PEW-3 

0.92 

0.010 

3.60 

1.47 

2.1 

PEW-4 

0 

0.001 

5.0 

0.20 

0.49 

Sum* 

11 

PEW-total" 

0.90 

0.013 

3.66 

1.41 

11 

Sum  represents  adding  mass  volume  from  individual  best  fit  for  each  extraction  well; 

‘  PEW-total  used  flux  average  mass  discharge  data  from  the  four  extraction  wells  to  fit  EST  model  and 
obtained  mass  volume  directly. 
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Table  5.  PLM  parameters  and  volume  of  mass  in  the  swept  volume  for  each  individual  extraction  well  and 
flux  averaged  of  all  wells.  The  results  are  the  best  fit  of  the  enhanced  dissolution  data  (mass 
discharge)  with  the  PLM  using  a  generic  algorithm.  Well  PEW-4  was  manually  fitted  because  the 
_ trend  of  data  is  not  obvious  and  best  fit  could  not  obtained  by  using  the  generic  algorithm. 


E 

/. 

r 

VW(L) 

PEW-l 

0.94 

0.016 

1.064 

1.9 

PEW-2 

0.81 

0.026 

0.94 

3.1 

PEW-3 

0.92 

0.010 

2.44 

2.2 

PEW-4 

0.0 

0.002 

0.01 

0.49 

Sum* 

7.7 

PEW-total” 

0.90 

0.013 

1.28 

7.8 

Sum  represents  adding  mass  volume  from  individual  best  fit  for  each  extraction  well; 

‘  PEW-total  used  flux  average  mass  discharge  data  from  the  four  extraction  wells  to  fit  PLM  model  and 
obtained  mass  volume  directly. 
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Table  6.  EST  model  parameters  used  to  predict  aqueous  dissolution 


fc 

&lnx 

^ampl(L) 

l^lnr 

0.010 

1.4 

2.0 

2.43 

0.010 

1.4 

20 

4.29 

0.010 

1.4 

11 

3.73 
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©  Injection  or  Extraction  Well 
©  Soil  Core 
©  MIP 


©  Mnlti-Lexel  Sampler 
☆  Pnsh-pnll  Tracer  Test  Location 


Figure  7.1.  DNAPL  source  zone  and  field  activities.  Soil  cores  at  locations  of  PEW-03,  SMLS-1,  SMLS- 
2  and  SIW-02  were  also  collected  at  locations  where  there  were  high  PID  signals.  Note:  dark  blue 
square  indicates  the  suspected  source  zone  area,  and  yellow  circle  shows  the  minimum  extent  of 
TCE  distribution  with  radius  1.2  m  and  thickness  1.4  m 
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Photoionization  detector  (PID)  signals  TCE  Concentration  (mg  TC'E/kg  Soil) 
(Max  Response  mV) 


Figure  7.2.  High  DNAPL  zone  (a)  MIP  -7;  (b)  Soil  core  data  at  location  SPW3-2;  (c)  pressure  from  a 
hydraulic  profiling  tool  (FIPT)  located  in  the  source  zone.  These  three  locations  are  close  and 
potentially  show  high  DNAPL  zone  at  the  depth  of  5.2  m  to  6. 1  m  within  a  low-flow  zone. 
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Figure  7.4.  Pore  velocity  (m/day)  versus  depth  for  MLSs  in  the  source  zone  during  the  PTT 
period.  Pore  velocity  was  calculated  based  on  the  distance  from  tracer  release  plane  (SIWs)  and 
nonpartitioning  tracer  (methanol)  peak  arrival  times  without  considering  tortuosity. 
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Figure  7.5.  Partitioning  tracer  breakthrough  curves  for  n-hexanol  and  2-octonal  with  exponential 
extrapolations  for  the  flow-weighted  average  PEW-T. 
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Figure  7.6.  Aqueous  TCE  concentration  during  the  enhanced  dissolution  and  after  the  enhanced 
dissolution  at  SMLS-1,  4,  7,  8  transect. 
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0  10  20  30  40  50 

Time  (days) 

Figure  7.7.  The  capture  time  estimation  for  flux  average  TCE  concentration  from  the  four 
extraction  wells.  The  dashed  line  represents  the  capture  time  (2.64days)  when  the  flow 
approached  steady  state  flow. 
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Figure  7.8.  EST  model  best  fit  for  enhanced  dissolution  data.  Note:  the  pore  volume  zero  represents  the 
date  for  groundwater  sampling  before  the  enhanced  dissolution. 
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Pore  Volume 


EST  prediction  (20  L) 

EST  prediction  (2  L) 

EST  prediction  (11L) 

A  data  from  PEWs  during  enhanced  dissolution 
O  PFM  data  after  enhanced  dissolution 


—  —  —  PLM  prediction  (20  L) 

—  •  •  PLM  prediction  (2  L) 

. PLM  prediction  (TIL) 

□  PFM  data  after  enhanced  dissolution 


Figure  7.9.  Aqueous  dissolution  predictions  by  combining  soil  core  data,  MIP  data,  Push-Pull 
tracer  data  and  PTT  data  for  both  EST  and  PLM  models.  Passive  flux  meter  (PFM)  data  before 
enhanced  dissolution  (diamond  symbol)  and  after  enhanced  dissolution  (square  symbol)  were 
plotted. 
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9.  Calf  Pasture  Point  Assessment  using  Spatial  and  Temporal 
Transects  within  the  Plume 


a.  Evaluation  of  Source  Strength  Function 

An  understanding  of  the  source- strength  function  is  important  because  of  its  use  in  predicting 
responses  to  various  remedial  approaches,  and  therefore  is  a  valuable  tool  for  site  management. 
In  the  Calf  Pasture  Point  site  study,  the  initial  contaminant  is  estimated  to  be  about  7,500  gallons 
of  decontaminating  agent  non-corrosive  (DANC)  buried  in  drums.  The  composition  of  DANC  is 
predominantly  1,1,2,2-tetrachloroethane  (PCA)  with  an  estimated  release  mass  of  45,000  kg,  but 
this  estimate  has  a  high  uncertainty.  PCA  typically  degrades  rapidly  to  trichloroethylene  (TCE), 
which  is  the  predominant  contaminant. 

Wells  are  located  throughout  the  site,  and  the  focus  is  on  the  wells  screened  across  the  lower  silt 
and  till  units  (D-series  wells)  since  they  contain  the  highest  concentrations.  An  analysis  of  the 
groundwater  flow  field  was  also  conducted  and  shown  in  Figures  8.1  and  8.2. 
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Figure  8.1.  Groundwater  elevation  contour  map  based  on  average  groundwater  elevation 
measurements.  Blue  lines  represent  contour  elevations  in  feet  above  mean  sea  level,  and  the 
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grey  arrows  represent  groundwater  flow  directions.  The  bold  red  line  represents  the  curvilinear 
SZCP. 


Figure  8.2.  Coefficient  of  variation  for  groundwater  elevation  measurements. 

The  source-zone  control  plane  (SZCP)  is  a  control  plane  that  measures  the  mass  discharge  of 
DNAPL  immediately  down  gradient  from  the  source  zone.  The  estimation  of  mass  discharge  as  a 
function  of  time  can  also  be  called  the  source- strength  function. 

The  location  of  the  SZCP  impacts  the  source-strength  function,  and  the  SZCP  has  high 
uncertainty  since  the  limits  of  the  source  zone  have  not  been  determined  through  source-zone 
characterization.  The  location  of  the  SZCP  was  based  on  wells  with  long-term  monitoring  data, 
the  proximity  of  the  wells  to  the  inferred  DNAPL  release  location,  and  the  location  of  the  wells 
relative  to  the  groundwater  flow  field.  The  spacing  of  the  wells  is  roughly  90  feet.  Well  spacing 
has  a  large  impact  on  the  uncertainty  of  the  estimated  source  strength  function.  The  larger  the 
well  spacing,  the  greater  the  uncertainty  will  be.  Variability  in  concentration  with  depth,  and 
groundwater  flux  with  depth  impact  the  contaminant  flux  distribution.  The  travel  time  was 
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estimated  to  be  20  years,  but  other  evidence  suggests  the  travel  time  is  most  likely  less  than 
estimated. 

Concentration-time  series  information  from  wells  04D,  05D,  and  17D  was  used  to  estimate  the 
source-strength  function.  Combining  this  information  with  well  specific  groundwater  flux, 
contaminant  flux  was  estimated  which  was  then  applied  to  calculate  the  mass  discharge.  The 
PCA  was  converted  to  TCE  on  a  molar  basis  to  calculate  mass  balance  since  TCE  is  the 
predominant  contaminant.  Mass  discharge  ranged  from  300  g/day  in  1996  to  50  g/day  in  2010. 
The  Power-Law  Model  was  used  to  related  DNAPL  mass  in  the  source  zone  to  mass  discharge 
from  the  source.  The  first  measurements  used  for  the  model  were  from  1996.  Certain  calculations 
from  October  2008  were  completed,  but  not  all  calculations  were  included. 

The  Power-Law  Model  can  also  be  used  to  estimate  mass  in  1968  by  extrapolating  data  from 
1996  under  the  assumption  that  the  model  does  not  change  with  time.  Two  cases  were 
considered:  the  parameter  T  =  0,  so  the  mass  discharge  over  42  years  would  be  4,600  kg,  and  the 
parameter  T  =  1,  so  the  mass  discharge  over  42  years  would  be  22,400  kg.  The  first  case  resulted 
in  a  mass  discharge  that  should  generally  be  three  times  greater.  The  second  case  is  more  likely 
within  the  plume  mass  estimates,  but  the  mass  discharge  in  1968  would  have  been  6,300  g/day 
which  is  higher  than  expected,  but  not  impossible. 

Mass  balance  requires  that  the  mass  discharged  from  the  source  zone  equal  the  mass  within  the 
plume,  plus  any  mass  that  may  have  been  degraded.  Plume  mass  estimates  are  useful  for  source 
strength  function  evaluation.  Sorbed  mass  is  accounted  for  using  an  estimate  of  the  site  specific 
retardation  factor,  R  =  7.  Since  PCA  degrades  to  TCE,  which  then  degrades  to  DCE,  the 
degradation  is  accounted  for  by  using  the  molar  sum  of  all  three  contaminants.  The  estimates  for 
the  source  zone  mass  are  roughly  100-300  kg  PCA,  and  the  plume  mass  is  roughly  13,000- 
20,000  kg  PCA,  TCE,  and  DCE.  The  ratio  of  contaminants  in  the  source  zone  to  the  mass  in  the 
plume  is  the  site  age,  and  thus  this  site  is  a  relatively  aged  site. 

Three  transects  were  formed  to  conduct  mass  flux  and  mass  discharge  calculations  based  on 
PFM  deployments.  Approximately  90  PFMs  were  deployed  in  three  transects  to  estimate  mass 
discharge.  The  contour  plots  of  mass  flux  are  provided.  These  were  used  to  calculate  mass 
discharge  at  each  transect  and  it  can  be  observed  that  mass  discharge  increases  along  the  length 
of  the  plume  thus  suggesting  an  exponential  decay  model  may  be  appropriate  for  the  source 
strength  function.  This  can  be  estimated  based  on  the  observed  data. 

The  longevity  of  the  DNAPL  estimated  to  remain  within  the  source  zone  can  be  evaluated  using 
a  first  order  degradation  analogy.  The  Power-Law  Model  using  T  =  1  is  a  good  estimate,  but  this 
does  not  include  other  factors  that  may  prolong  contaminant  mass  in  the  source  zone. 

Recommendations  for  future  work  are  also  included  in  the  next  section.  The  recommendations 
are  designed  to  minimize  uncertainty  with  this  analysis,  and  future  site  management  decisions. 
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Figure  8.3.  Source  strength  functions.  Measured  data  shown  in  blue  diamonds,  model 
fits  shown  for  T  =  {0.5,  1.0,  2.0}  as  green,  purple,  and  brown  lines. 


Calf  Pasture  Point  -  Installation  of  Control  Planes 


Figure  8.4.  Three  transects  at  Calf  Pasture  Point  for  PFM  deployments. 
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Control  Plane  at  MW07-39D 
Deep  and  Intermediate  Depths 
Equivalent  TCE  Flux  (g/sqm/day) 


SOURCE 


Figure  8.5.  Three  transect  mass  flux  based  on  PFM  deployments. 


Figure  8.6.  Source  strength  function  based  on  three  transects  sampled  at  Calf  Pasture 
Point. 
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b.  Uncertainty  Implications  for  Predicting  Source  Dissolution 

The  contaminated  site  is  assumed  to  be  primarily  PC  A  with  an  estimated  release  mass  of  45,000 
kg.  Due  to  the  assumptions  and  estimate,  the  release  mass  is  expected  to  have  high  uncertainty. 
The  source  strength  function  uses  historical  data  and  field  measurements  to  estimate  mass 
discharge  as  a  function  of  time.  The  source  strength  function  is  a  model  of  the  data  used  to 
predict  responses  to  remedial  approaches,  but  it  can  contain  uncertainty. 

Since  the  limits  of  the  source  zone  have  not  been  determined,  the  source  zone  area  is  subject  to 
uncertainty.  Some  wells  were  also  not  included  in  the  long-term  monitoring  program.  The  wells 
that  were  selected  were  five  wells  that  were  about  450  feet  in  distance,  and  they  had  a  well 
spacing  of  90  feet.  The  spacing  of  the  wells  impact  the  accuracy  of  the  source  strength  function, 
and  uncertainty  arises  from  the  large  unsampled  space  between  the  wells.  The  vertical 
distribution  of  the  contaminants  near  the  source  zone  was  not  evaluated,  but  is  important  for  the 
source  zone  control  plane.  Therefore  the  vertical  depth  was  assumed  to  equal  the  screened 
interval  of  the  wells. 

The  Power-Law  Model  was  used  to  predict  the  mass  in  1968.  It  was  assumed  the  model  did  not 
change  with  time;  however,  that  assumption  is  most  likely  not  correct,  and  contains  uncertainty. 

Plume  analysis  is  useful  to  estimate  mass  balance  between  the  plume  and  the  source  zone.  Mass 
discharge  from  the  source  zone  must  equal  the  mass  in  the  plume,  as  well  as  any  degraded  mass. 
Sorbed  mass  can  be  estimated  using  the  site  specific  retardation  factor,  R.  However,  the 
uncertainty  in  the  plume  mass  estimate  also  becomes  a  function  of  the  uncertainty  in  R. 
Typically,  plume  mass  estimates  have  a  high  uncertainty  which  is  most  likely  from  the  sparse 
sampling  density. 

DNAPL  source  longevity  can  be  evaluated  using  a  first-order  degradation  analogy.  These 
estimates  do  not  include  factors  that  may  prolong  contaminant  mass  in  the  source  zone,  such  as 
sorption  to  the  solid  matrix  and  diffusion  from  low  conductivity  zones,  and  thus  they  have  a  high 
uncertainty. 

Recommendations  for  minimizing  the  uncertainty  in  analyzing  source  dissolution  and  assessing 
site  management  decisions  are  provided.  Wells  within  the  source  zone  area  that  are  not  currently 
part  of  the  long-term  monitoring  program  should  be  sampled  in  order  to  better  evaluate  the 
strength  and  distribution  of  the  source.  It  is  also  recommended  that  a  preliminary  passive  flux 
meter  deployment  be  made  to  evaluate  the  vertical  distribution  of  the  contaminants.  This  would 
also  provide  information  to  improve  the  source  strength  function.  Following  the  passive  flux 
meter  deployment,  additional  soil  borings  are  to  be  installed  within  the  source  zone.  They  will  be 
analyzed  for  chlorinated  ethane  and  ethane  compounds,  as  well  as  organic  carbon.  After 
additional  wells  are  installed,  it  is  proposed  that  a  second  passive  flux  meter  deployment  be 
conducted. 


10.  Back-Diffusion  at  Field  Sites  Using  Flux  Measurements 

The  former  Brandywine  Defense  Reutilization  and  Marketing  Office  (DRMO)  Yard  in 
Brandywine,  MD,  is  and  inactive  federal  facility  that  was  in  operation  from  1943-1987.  It  was 
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used  for  temporary  storage  of  scrap  materials  and  hazardous  waste,  and  is  now  on  the  National 
Priorities  List  (NPL).  The  DRMO  yard  groundwater  is  contaminated  primarily  with  TCE,  and 
the  plume  which  is  underlain  by  the  Calvert  Foundation  aquitard  has  extended  past  the  property 
into  a  residential  neighborhood.  In  2006  an  Interim  Record  of  Decision  (ROD)  was  signed,  and 
in  2008  the  ROD  was  implemented  to  capture  the  3-acre  core  of  the  plume.  The  remaining  20- 
acre  TCE  plume  was  treated  with  In-situ  Bioremediation  (ISBR).  38.6  kg  of  VOCs  were 
removed  during  5  years  of  treatment. 
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Figure  9.1:  Site  map  showing  pre-remedy  extent  of  TCE  plume  (upper  left),  ISBR  treatment 
areas,  and  post-remedy  implementation  TCE  plume  extent  (March  2013). 

The  Brandywine  Formation  and  the  Calvert  Formation  are  located  beneath  the  DRMO  yard.  The 
upper  Calvert  is  21-30  feet  below  ground  surface,  and  30-42  feet  below  ground  surface  is  the 
oxidize  Calvert.  Both  are  considered  the  aquitard.  The  Brandywine  Formation  contains  the 
contaminant  plume  in  the  upper  12  feet,  and  is  comprised  of  clay,  silt,  sand,  and  gravel.  The 
Brandywine  Formation  has  been  determined  to  be  heterogeneous  laterally  and  vertically  across 
the  site,  and  consists  of  four  distinct  strata:  Shallow  Brandywine,  Upper  Intermediate 
Brandywine,  Fower  Intermediate  Brandywine,  and  Deep  Brandywine. 

Investigation  suggested  that  limited  source  material  remained  within  the  presumed  release  area; 
however  a  study  revealed  contaminant  source  remained  when  the  pump  and  treat  system  were 
shut  down.  Membrane  interface  probe  technology  was  selected  to  collect  real-time,  continuous 
data  on  the  distribution  of  volatile  organic  compounds.  The  area  selected  for  investigation  was 
based  on  VOC  results  for  soil  samples.  Combined  with  the  MIP,  there  is  a  conductivity  logging 
tool  that  interprets  lithology.  The  MIP  is  used  to  identify  zones  of  contamination.  There  were 
three  electric  conductivity  logging  tools  tested:  the  flame  ionization  detector  (FID),  the 
photoionization  detector  (PID),  and  the  electron  capture  detector  (ECD).  The  ECD  was  chosen 
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for  analysis  of  the  site  since  it  provides  the  lowest  possible  detection  levels  for  chlorinated 
compounds. 


After  the  MIP  investigation,  5  monitoring  wells  with  15  foot  screens  were  installed.  Flux  well 
(FW)-l  was  installed  in  the  presumed  release  area,  FW-3  was  installed  down-gradient,  FW-5 
was  installed  further  down-gradient  where  it  was  presumed  contamination  was  not  prevalent,  and 
FW-2  and  FW-4  were  installed  side-gradient  which  created  a  transect  with  FW-3.  Passive  Flux 
Meters  (PFM)  were  deployed  in  each  well  to  determine  the  mass  flux  of  water  and  contaminant 
flowing  per  unit  area  in  a  given  period  of  time.  Three  5-foot  PFMs  were  stacked  across  the  entire 
screen  length  for  20  days.  Once  they  were  retrieved  they  were  sampled  at  1-foot  intervals. 
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Figure  9.2:  Site  map  showing  MIP  boring  and  flux  well  locations. 


The  ECD  results  indicated  that  the  high  concentrations  of  TCE  are  located  in  the  transitional  silt 
layer,  and  the  highest  1  meter  into  the  oxidized  Calvert.  The  results  also  revealed  that  down- 
gradient  from  the  release  area,  and  west  of  the  groundwater  extraction  trench,  the  contamination 
is  within  the  groundwater,  and  not  the  transitional  silt  layer. 
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Electrical  C  onductivity 


Figure  9.3:  ECD  and  conductivity  MIP  responses  along  the  axis  of  the  plume  A-A\ 
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Figure  9.4:  Lateral  extent  of  TCE  contamination  determined  by  the  ECD. 


The  TCE  flux  profiles,  determined  from  the  PFM  results,  show  that  the  5  wells  contain 
significant  variability  and  support  the  conceptual  model  of  high  mass  flux  associated  with  the 
aquitard.  FW-3  contained  the  highest  mass  flux  and  the  most  significant  contrast  in  the  center  of 
the  suspected  back-diffusion  zone  within  the  plume  core. 
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Figure  9.5:  Flux  well  profiles  showing  Darcy  velocity  and  TCE/DCE  flux. 


The  highest  mass  flux  was  located  at  a  depth  corresponding  to  the  transitional  silt  interface.  The 
two  side-gradient  wells  contained  lower  mass  flux.  The  MIP  profile  for  FW-2  reveals  a  distinct 
mass  source  22-35  feet  below  ground  surface  that  is  likely  providing  mass  flux  to  the 
Brandywine  Formation.  There  was  also  a  slight  increase  in  the  Darcy  velocity  in  the  flux  profile 
that  correlates  with  a  lower  conductivity  response.  In  FW-2,  the  highest  concentration  observed 
was  located  at  the  upper  screened  portion  of  the  well  which  corresponds  to  the  top  of  the  TCE 
contaminated  zone.  The  result  of  this  supports  the  conceptual  model  of  mass  diffusing  from  the 
silt  interface  into  the  Brandywine  Formation. 
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Figure  9.6:  FW-2  flux  profile  compared  to  MIP  ECD  and  conductivity  responses. 
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FW-5  results  in  a  much  lower  mass  flux  that  FW-3,  and  also  reveals  a  shift  from  TCE  to  DCE. 
The  MIP  profile  shows  low  mass  stored  in  the  transitional  silt,  and  suggests  that  mass  may  still 
be  loading  the  aquitard. 


Table  1.  Mass  discharge  per  unit  width  of  aquifer  for  each  well. 


Darcy 

TCE 

1,2-DCE 

Well 

(cm/day) 

(mg/m/day) 

(mg/m/day) 

FW-1 

1.3 

73.8 

3.8 

FW-2 

1.5 

47.5 

3.3 

FW-3 

1.1 

891.1 

71.8 

FW-4 

1.2 

12 

8.5 

FW-5 

1.4 

15.2 

79.1 

The  mass  discharge  was  calculated  using  the  summation  of  the  product  of  flux  and  the  cross- 
sectional  flow  area.  The  PFM  based  mass  flux  results  combined  with  the  MIP  ECD  profile  in  FW- 
3  provide  the  ability  to  link  mass  flux  from  back-diffusion  with  mass  flux  observed  in  the  aquifer. 
A  simple  approach  to  making  this  link  is  to  use  a  simplified  one-dimensional  analytical  solution 
for  back  diffusion.  The  model  is  provided  here  and  assumes  a  step  input  of  constant  concentration 
to  the  diffusive  zone  over  a  period  of  time  specified  as  t*. 


r  \ 

L 1 

This  equation  provides  the  diffusive  flux  from  the  aquitard,  Jbd,  at  some  time  of  interest  t  given 
an  estimate  of  the  diffusion  coefficient,  D,  and  the  applied  loading  concentration  Co.  A  quick 
estimate  was  made  assuming  a  Co  value  of  100  mg/L  (10%  of  TCE  solubility)  a  diffusion 
coefficient  for  TCE  in  water  (1x10  9  m2/s)  and  a  loading  period  t*  of  20  years  and  5  years  after 
back-diffusion  started  (t=25  years).  The  calculated  flux  from  the  aquitard  is  0.004  g/m2/day.  If 
this  value  is  applied  over  a  length  of  50  m  the  mass  discharge  per  unit  width  of  aquifer  is  about 
0.35  g/m/day  which  is  of  the  same  magnitude  as  observed  in  well  FW-3  (0.89  g/m/day).  This 
supports  the  back-diffusion  conceptual  model  for  the  site. 

The  results  of  this  study  were  used  to  characterize  the  nature  and  extent  of  source  area 
contamination.  The  MIP  results  and  PFM  results  correlate  with  respect  to  contaminant  mass  within 
the  transitional  silt,  the  locations  and  magnitude  of  diffusive  mass  flux,  and  Darcy  velocity.  The 
characteristics  of  the  site  allow  the  plume  to  persist  in  the  Brandywine  Formation,  and  the  result 
of  the  study  has  provided  focus  on  the  evaluation  of  remedial  activities  and  which  areas  to  target. 
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11.  Conclusions  and  Implication  for  Future  Research 


This  project  focused  on  evaluating  DNAPL  source  strength  functions  based  on  field  data 
observations  and  the  use  of  simple  models.  In  this  effort,  several  methods  were  proposed  and 
evaluated  as  approaches  that  can  be  used  to  estimate  parameters  for  DNAPL  source  strength 
functions.  The  results  of  the  project  have  been  disseminated  through  a  series  of  journal  articles 
that  are  listed  in  Appendix  B.  The  research  started  with  formulation  of  source  strength  functions 
and  appropriate  models  and  parametrization  (Chen  and  Jawitz,  2009;  Wang  et  al.,  2013).  The 
work  then  was  applied  to  field  sites  with  adequate  data  allowing  determination  of  model 
parameters  to  characterize  source  strength  functions  for  DNAPL  sites  (Johnston  et  al.,  2013, 
2014;  Wang  et  al.,  2014).  The  evaluation  of  uncertainty  associated  with  determining  mass 
discharge  and  source  strength  function  was  evaluated  for  both  pumping  based  measurements 
(Xiaosong  et  al.,  2014)  and  point  level  measurements  (Brooks  et  al.,  2015).  Thus  guidance  in 
data  needs  and  approaches  were  described  for  those  interested  in  understating  DNAPL  source 
zone  function.  On  a  parallel  effort,  research  was  conducted  to  better  understand  the  fundamentals 
of  back  diffusion  from  low  permeability  media  in  the  subsurface.  This  started  with  a  theoretical 
approach  linking  DNAPL  source  strength  functions  to  back  diffusion  (Brown  et  al.,  2013).  Then 
a  series  of  laboratory  scale  studies  were  conducted  developing  light  transmission  techniques  for 
clays  (Yang  et  al.,  2014)  and  the  effects  of  different  aquifer  aquitard  layers  (Yang  et  al.,  2015) 
and  finally  linking  to  source  strength  functions  (Yang  et  al.,  2016).  Back  diffusion  at  field  scales 
was  then  assessed  (Yang  et  al.,  2017).  The  project  final  report  contains  some  of  the  material 
presented  in  these  publications  but  also  some  that  is  currently  being  revised  for  publication  such 
as  the  Calf  Pasture  Point  study.  The  final  report  and  associated  journal  articles  will  be  included 
as  part  of  an  EPA  Guidance  Document  of  DNAPL  source  strength  functions. 

While  some  links  between  plume  conditions  and  DNAPL  source  strength  functions  have  been 
evaluated  as  part  of  this  work,  this  is  an  area  that  could  be  greatly  improved.  Linking  conditions 
in  the  plume  to  the  complex  behavior  of  DNAPL  source  zones  is  challenging  given  that  most 
source  zones  have  a  history  of  partial  remedial  efforts.  This  complicates  our  understanding  of 
the  nature  of  the  plume  distribution.  While  Calf  Pasture  Point  was  ideal  due  to  the  lack  of 
remedial  efforts,  this  is  the  exception  and  better  understanding  of  more  complex  source  zone 
behavior  is  needed.  This  improved  understanding  of  DNAPL  source  zone  behavior  and  plume 
response  can  be  then  connected  to  the  assimilative  capacity  of  the  aquifer  and  alternative  mass 
discharge  or  mass  flux  objectives  can  be  formulated. 
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Appendix  A:  Supporting  Data 

Table  A1 .  Summary  of  mass  discharge  measurements  (g/day)  using  Method  1  at  Calf  Pasture  Point. 
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Table  A2.  Summary  of  mass  discharge  measurements  (g/day)  using  Methods  2a  and  2b. 
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Table  A3.  Summary  of  mass  discharge  measurements  (g/day)  using  Method  3. 
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Summary  of  Mass  Discharges  for  PFM3  and  PFM4 
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